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ABSTRACT 


An  experimental  technique  is  developed  for  measuring 
absolute  neutron  detection  efficiencies  and  subsequently 
used  in  the  structure  study  of  Na2 1  via  the  Ne2 0  (d,n) Na2 1 
reaction. 

Angular  distributions  of  the  neutrons  leading  to  the 
four  bound  states  of  Na21  were  obtained  at  a  deuteron  bom¬ 
barding  energy  of  3.924  MeV.  A  D.W.B.A.  analysis  was  used 
to  check  the  spin  and  parity  assignments  and  to  obtain  abso¬ 
lute  spectroscopic  factors  for  the  Na21  0.335  (5/2+) ,  1.751 
(7/2+) ,  and  2.433  (l/2+)  levels. 

A  strong  coupling  collective  model  interpretation  of 
the  Na21  level  structure  is  given  and  a  band-limited  sum¬ 
mation  rule  is  used  to  estimate  the  compound  reaction  contri¬ 
bution  to  the  formation  of  the  0.335,  1.751,  and  2.433  MeV. 


levels . 
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CHAPTER  I 

INTRODUCTION 

In  low  energy  nuclear  spectroscopy  one  has  two  basic 
interests.  One  is  to  elucidate  the  structural  content  of 
the  nucleus  of  interest,  and  the  other  is  to  determine  the 
nature  of  the  reaction  mechanism.  Progress  in  both  these 
areas  is  governed  by  the  extent  to  which  a  comparison  can 
be  made  between  one's  experimental  findings  and  various 
model  predictions. 

A  number  of  important  nuclear  parameters,  such  as  spin, 
parity,  and  excitation  energy  can  generally  be  determined 
on  the  basis  of  relative  cross  section  data.  However,  these 
parameters  tend  to  represent  the  more  gross  features  of  the 
nucleus  and  provide  only  limited  insight  into  its  internal 
make  up.  As  one  attempts  to  probe  more  deeply  into,  not 
only  the  internal  makeup  of  the  nucleus  but  also  into  the 
nature  of  the  reaction  mechanism,  a  more  quantitative  analysis 
is  required.  That  is,  only  through  absolute  cross  section 
data  can  one  attempt  to  answer  questions  such  as: 

(i)  to  what  extent  is  the  residual  nucleus 

represented  by  the  target  nucleus  plus  the 
captured  particle; 
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(ii)  to  what  extent  can  the  nucleon-nucleon 
interactions  be  represented  by  averaged 
potentials;  and 

(iii)  to  what  extent  is  there  competition  between 
the  direct  reaction  and  compound  reaction 
mechanisms? 

All  three  questions  relate  directly  to  the  two  fundamental 
problems  in  nuclear  physics:  that  of  determining  the  nature  of 
the  nuclear  forces,  and  that  of  treating  the  many  body  problem. 

In  this  thesis  a  report  is  given  of  work  undertaken  at  the 
University  of  Alberta  to:  (i)  develop  an  experimental  technique 
for  measuring  absolute  neutron  detection  efficiencies  and  (ii) 
apply  the  results  of  such  a  technique  to  the  absolute  cross- 
section,  structure  study  of  Na2  1  via  the  Ne 2 0  (d  ,n) Na 2 1  reaction. 

For  the  absolute  efficiency  measurement  the  D(d,n)He3 
reaction  was  used,  with  the  He3  particles  being  detected  in 
coincidence  with  their  associated  neutrons.  The  coincidence 
requirements  were  such  that  a  neutron  was  detected  if  and  only 
if  an  associated  He3  particle  was  also  detected.  However,  all 
He3  particles  were  detected  independent  of  their  associated 
neutrons.  The  ratio  between  the  number  of  neutrons  detected 
to  the  number  of  He3  particles  is  then  a  measure  of  the  absolute 
neutron  detection  efficiency. 

The  study  of  the  Ne2 0  (d ,n) Na2 1  reaction  was  initiated  as 
a  search  for  new  low  lying  structure  within  the  Na2 1  nucleus 
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(viz.  a  questionable  1.42  MeV  level  (Ma  56)  and  possible 
doublet  structure  in  the  0.335  (Gi  65)  and  2.433  (Po  67)  levels) 
and  to  obtain  absolute  cross  section  data  for  the  formation  of 
the  four  bound  states  of  Na21  (i.e.:  ground  (3/2+),  0.335  MeV 
(5/2+) ,  1.751  MeV  (7/2+),  and  2.433  MeV  (l/2+)  levels).  The 
Ne20  target  consisted  of  a  nickel-windowed  gas  cell  and  the 
neutrons  were  detected  via  a  time-of-f light  technique.  Both 
structure  and  reaction  interpretations  of  the  Na2 1  nucleus  are 
made  within  the  framework  of  the  strong  coupling  collective 
model . 
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CHAPTER  II 

EXPERIMENTAL  DETAILS 

2.1  Introduction 

As  a  consequence  of  the  neutron's  neutral  charge  there 
is  no  unambiguous  relationship  between  the  recoil  pulse  height 
spectrum  of  a  neutron  scintillator  and  the  energy  of  the 
incident  neutron.  Thus,  to  date,  the  most  precise  method  of 
measuring  neutron  energies  is  to  determine  particle  velocities 
by  individually  timing  their  flight  over  an  accurately  known 
distance.  Such  a  technique  requires  the  accurate  determi¬ 
nation  of  two  time  points  corresponding  to  the  generation  and 
detection  of  the  neutron.  These  time  points  are  most  accu¬ 
rately  determined  by  using  a  pulsed  beam  in  conjunction  with 
fast  timing  circuits.  Such  a  technology,  in  recent  years,  has 
developed  into  a  very  powerful  tool  available  for  fast  neutron 
spectroscopy  in  the  nano-second  time  scale.  The  significant 
features  are: 

(i)  good  energy  resolution 

(ii)  high  detection  efficiency 

(iii)  applicability  to  a  wide  range  of  neutron  energies. 
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The  start  signal  originates  from  a  pick  off  capacitor 
located  within  the  beam  tube  next  to  the  target.  The  beam 
pulse  induces  a  charge  which  is  amplified  and  shaped 
appropriately . 

The  corresponding  stop  signal  originates  from  the  neutron 
detector,  which  consists  of  an  organic  scintillator  optically 
coupled  to  a  photomultiplier  tube. 

The  start-stop  signals  are  then  fed  into  a  fast  time-to- 
amplitude  converter  and  the  corresponding  output  signal  is 
analyzed  by  a  multichannel  pulse  height  analyzer. 

The  length  of  the  beam  pulses  are  short  in  comparison  to 
the  corresponding  neutron  flight  times.  This  serves  to  maxi¬ 
mize  energy  resolution  which,  for  non-relativistic  neutrons, 
is  given  by 


(where  At  is  the  width  of  the  pulse,  t  the  flight  time  of  the 
neutron,  E  the  energy  of  the  neutron,  and  AE  the  resulting 
uncertainty  in  the  neutron  energy.)  In  addition  the  pulse 
repetition  rate  is  sufficiently  slow  to  allow  the  low  energy 
neutrons  to  be  detected  before  the  arrival  of  the  next  pulse. 

For  the  system  employed  at  the  University  of  Alberta  an 
intrinsic  pulse  width  of  less  than  0.5  ns.  is  obtainable  at  a 
repetition  rate  of  1  mH.  The  neutron  flight  path  is  variable 
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from  1  to  6  meters.  And  the  multichannel  pulse  height  analyzer 
consists ,  in  part,  of  an  on  line  SDS-920  computer. 

2 . 2  Beam  Transport  System 

Fig.  2-1  is  a  schematic  representation  of  the  beam  trans¬ 
port  system  used  at  the  University  of  Alberta  for  the  generation 
of  a  high  intensity  pulsed  beam.  The  system  employs  both  pre¬ 
acceleration  pulsing  and  post-acceleration  compression,  a  combi¬ 
nation  that  serves  to  maximize  the  beam  current  without  increasing 
the  pulse  width. 

A  d.c.  beam  of  several  milliamps  is  generated,  within  a 
source  bottle,  at  the  high  voltage  terminal  of  the  Van  de 
Graaff^  by  a  125  mH.  source  oscillator  in  conjunction  with  a 
source  magnet.  The  beam  is  brought  out  through  an  opening 
known  as  the  canal  and  focused,  by  a  set  of  Einzel  lens  and  a 
terminal  analyzing  magnet,  at  a  1/16  inch  aperture.  Pre¬ 
acceleration  pulsing  is  achieved  by  sweeping  the  beam,  in  an 
elliptical  pattern,  over  the  1/16  inch  aperture  by  a  pair  of 
R.F.  deflection  plates  set  at  right  angles  to  each  other  and 
with  their  voltages  90°  out  of  phase.  The  R.F.  frequency  is  set 
at  1  mH.  producing  a  triangular  pulse  of  10  to  15  ns.  pulse 
width.  Additional  focusing  is  achieved  by  another  lens, 
the  self  focusing  column,  and  a  quadruple  doublet. 

^  A  6.00  MeV  Van  de  Graaff  Accelerator,  Model  CN,  High 

Voltage  Engineering  Corp. ,  Burlington,  Massachusetts,  U.S.A. 
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FIGURE  2-1  BLOCK  DIAGRAM  OF  BEAM  TRANSPORT  SYSTEM 
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At  the  base  of  the  accelerator  is  a  90°  mass-energy 
analyzing  magnet  which  focuses  the  beam  at  a  pair  of  energy 
stabilizing  slits  some  57  inches  away.  The  slits  are  in  turn 
servo^connected  to  a  corona  discharge  probe.  The  varying 
current  load  of  the  discharge  probe  serves  to  stabilize  the 
beam  energy. 

Next  a  magnetic  quadruple  lens  refocuses  the  beam,  after 
which  it  passes  between  the  plates  of  a  10  mH. ,  R.F.  deflection 
oscillator  which  sweeps  the  beam  across  the  entrance  aperture 
of  a  90°  Mobley  compression  magnet.  The  oscillator  is  phase 
locked  to  the  ion  burst  via  a  cylindrical  pick-off  probe  placed 
within  the  beam  tube  just  before  the  R.F.  unit.  Time  compression 
is  achieved  by  sweeping  the  beam  over  decreasing  path  lengths 
which  have  a  common  focal  point  at  the  target.  Since  the  leading 
particles  have  a  longer  flight  path  they  tend  to  arrive  at  the 
target  coincident  with  the  trailing  particles.  The  beam  bunching 
factor  of  the  present  installation  is  of  the  order  of  12 ,  giving 
an  intrinsic  pulse  width  of  less  than  0.5  ns. 

The  primary  factor  that  determines  the  bunching  performance 
of  the  magnet  is  known  as  the  beam  emittance.  It  is  a  measure 
of  the  transverse  dimensions  and  angular  divergence  of  the  pre¬ 
bunched  pulse.  Its  importance  and  relation  to  the  design 
criteria  of  the  Mobley  magnet  are  discussed  in  detail  in  the 
Ph.D.  thesis  of  W.  G.  Davies  (Da  66). 


'  ■ 

,  '  ■  1  .  1  1 

>  ■  •  -  -•  ;  ' 

•  J  "  -  -  •  ■  : 

:  .  :  -  '  Ni 

i<  nj* *•  )t  .3  »  J  .  3  stood  j  a.  edffct  n  • . .!  3tf7  '  ;  •  * 

1 

‘ 

‘ 

'  ■  .  ’  vi  :  Ci  HO d'  6  fTT 

:  : 

•  '.-‘I  ' 

*• 


-8- 


2 . 3  Time  of  Flight  Electronics 

Fig.  2.2  is  a  schematic  representation  of  the  electronics 
making  up  the  time-of-f light  spectrometer.  Parallel  to  the 
main  detector  is  a  monitor  system  from  which  a  normalization 
factor  is  obtained  for  angular  distribution  studies. 

The  main  detector  consists  of  a  Philips  XP  1040  photo¬ 
multiplier  tube  optically  coupled  to  a  cylindrical  (3.45"  x 

0.75")  quartz  container  filled  with  Ne  213^  liquid  scintil- 

2 ) 

lator.  The  monitor  detector  consists  of  a  Naton  phosphor 
optically  coupled  to  a  RCA  8575  photomultiplier  tube. 

The  main  detector  is  mounted  on  a  movable  cart  which, 
under  control  from  the  console,  can  be  positioned  in  both  detec¬ 
tion  angle  and  flight  path.  The  monitor  is  mounted  on  a  tripod 
and  generally  held  at  a  fixed  position. 

The  stop  signal,  for  both  systems,  comes  from  a  3  inch 
long  cylindrical  capacitor  placed  within  the  beam  tube  some  30 
cm.  in  front  of  the  target.  Before  reaching  the  stop  channel  of 
the  time-to-amplitude  converter  (T.A.C.)  the  pulse  is  delayed  suffi 
ciently  to  allow  a  start  pulse,  corresponding  to  the  lowest  energy 
neutrons,  to  reach  the  T.A.C.  first.  By  so  inverting  the  time 
sequence  one  can  decrease  the  overall  dead  time  of  the  system. 

This  is  a  consequence  of  the  virtual  elimination  of  those 

^Nuclear  Enterprises  Ltd.,  Winnipeg,  Manitoba. 

2 ^ Nash  &  Thompson,  Ltd.,  Hookrise  South,  Tolworth,  Surrey, 
England . 
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saturation  pulses  corresponding  to  a  start  signal  but  no  stop. 
That  is,  in  the  event  that  no  neutrons  are  detected  for  a  given 
ion  burst  the  spectrometer  is  not  activated  but  remains  ready 
for  the  next  beam  burst. 

The  time-to-amplitude  converters  are  based  on  a  design  by 
D.  L.  Wieber  (Wi  63)  using  capacitance  charging  and  a  time 
base  of  about  300  ns.  Cross  over  timing  (Ge  67)  is  used  in  both 
channels  of  the  monitor  system  but  only  in  the  stop  channel  of 
the  main  system.  For  the  main  detector  a  signal  is  taken  from 
the  14th  dynode  of  the  photo  tube,  inverted,  and  fed  into  a 
start  discriminator  shaper,  designed  by  Dr.  G.  C.  Neilson.  The 
leading  edge  of  this  pulse  serves  as  the  start  signal.  A 
second  pulse  taken  from  the  11th  dynode  of  the  photo  tube 
serves  as  a  linear  channel  from  which. 

(i)  a  discrimination  level  is  determined  via  signal 
channel  analyser  No.  3,  and 

(ii)  a  time  base  correction  is  made  via  ADC-B  and  the 
on  line  computer. 

The  single  channel  analyzer  is  operated  in  the  discrimi¬ 
nation  mode  only  and  the  bias  level  is  determined  as  a  compro¬ 
mise  between  neutron  threshold  energy  and  detection  efficiency 
on  the  one  hand  and  energy  resolution  and  background  on  the 
other.  A  threshold  energy  of  0.80  MeV.  was  used  throughout  this 
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The  digital  information  entering  the  computer  via  ADC-B 
is  a  measure  of  the  pulse  height  of  the  start  signal  activating 
T.A.C.  No.  2.  Since  time  zero  is  determined  by  when  this 
pulse  crosses  a  fixed  bias  level,  varying  time  points  will 
correspond  to  varying  pulse  heights.  However,  once  having  deter¬ 
mined  the  relationship  between  pulse  height  and  time  shift  the 
knowledge  entering  the  computer  via  ADC-B  allows  one  to  correct 
for  this  time  walk,  resulting  in  an  increase  in  resolution 
(0.9  ns.  from  1.25  ns.). 

Such  a  correction  is  achieved  by  a  general  purpose  program 
(GPKS)  written  by  Dr.  W.  K.  Dawson.  First  the  pulse  height 
spectrum  of  a  mono-energetic  neutron  peak  is  divided  into  64 
bins.  The  time  shift  for  each  bin  is  determined  and  entered  as 
a  table  in  the  computer.  Thus  every  pulse  that  enters  the  com¬ 
puter  via  ADC-B  is  associated  with  one  of  the  64  bins.  The 
corresponding  time  correction  is  determined  from  the  table  and 
added  to  the  timing  data  entering  the  computer  via  ADC-A. 

Both  analog  to  digital  counters^ ,  ADC-A  and  ADC-B  are 
1024  channel  units,  which  along  with  the  walk  correction  opera¬ 
tion  produces  a  system  dead  time  of  the  order  of  700  ys.  This 
allows  a  maximum  conversion  rate  of  about  1500  counts  per 
second.  During  the  time  GPKS  is  handling  data,  input  to  the 
computer  is  inhibited  by  the  coincidence  interrupt  unit.  Its 

^  Technical  Measurement  Corporation,  Model  217A,  North 
Haven,  Conn.,  U.  S.  A. 
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detailed  operation  is  discussed  in  an  internal  report  by  D.  A. 
Gedcke  (Ge  66) . 

Signals  from  the  anode  and  13th  dynode  of  the  photo  tube 
provide  inputs  to  a  y-ray,  neutron  discriminator  of  the  pulse- 
overlap  design  (Re  66).  The  output,  which  signifies  the 
identification  of  the  pulse  as  a  neutron  is  held  in  coincidence 
with  signals  validating  it  as  being  above  the  neutron  threshold 
level  and  also  initiating  a  legitimate  start  pulse.  The  output 
of  the  triple  coincidence  is  then  fed  into  a  routing  unit 
designed  by  D.  A.  Gedeke  (Ge  66a)  which  routes  the  corresponding 
corrected  timing  data  into  a  second  1024  channel  memory  range  of 
the  computer.  The  first  1024  channel  range  then  contains  those 
data  that  the  n-Y  discriminator  identified  as  y-rays. 

In  addition  to  storing,  correcting,  and  displaying  data 
the  utility  of  GPKS  allows  one  to  transfer  the  memory  contents 
of  the  T.M.C.  210  unit,  which  contains  the  associated  monitor 
spectrum,  to  the  memory  of  the  computer.  All  data  is  then 
easily  transferred  to  magnetic  tape  for  storage  purposes. 

2 . 4  Gas  Target 

Fig.  2-3  is  a  design  drawing  of  the  gas  target  cell  used 

20  21 

throughout  most  of  the  Ne  (d,n)Na  work. 

The  cell  proper  consists  of  a  sheet  of  0.010"  tantalum 
folded  into  a  rectangular  parallelepiped  of  dimensions 
l1/4  H.  x  1"  W.  x  1/4"  T. 


The  enclosure  is  capped,  top  and 
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FIG.  2-3:  GAS  TARGET  CELL 


-12- 


bottom,  with  two  slotted  copper  discs.  An  epoxy"'")  is  used 

to  seal  together  the  overlapping  sides  of  the  parallelepiped 

and  the  copper-tantalum  interface,  resulting  in  a  leak  tight 

enclosure.  A  3/8"  hole  is  then  made  in  one  face  of  the 

enclosure  and  covered  over  by  a  50  yin.  nickel  foil.  A  second 

2 ) 

grade  of  epoxy  is  used  to  seal  the  foil  in  place. 

The  whole  assembly  screws  into  an  air  cooled  heat  sink 
which  serves  to  dissipate  the  heat  produced  by  the  beam  being 
stopped  within  the  gas  cell.  The  heat  sink  is  electrically 
insulated  from  the  target  assembly  by  a  teflon  washer,  which 
allows  one  to  monitor  the  charge  incident  on  the  target. 

Gas  is  fed  to  the  cell,  under  controlled  pressure,  via  an 
external  gas  handling  system.  This  pressure  data,  along  with 
the  temperature  data  obtained  from  a  thermocouple  connected  to 
the  heat  sink,  allows  for  an  accurate  determination  of  the 
target  thickness. 

The  system  functioned  well  at  a  bombarding  energy  of  4.00 
MeV. ,  a  gas  pressure  of  45  cm.  of  Hg . ,  and  a  beam  current  of 
1.5  yamps  with  neither  the  nickel  foil  nor  the  tantalum  backing 
contributing  significantly  to  the  background.  However,  when 
attempts  were  made  to  raise  the  bombarding  energy  to  5.00  MeV. 

Plastic  Steel,  Devcon  Canada  Ltd. 

2 ) 

Epoxy  Glue,  Borden  Chemical  Co. 
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the  foils  continued  to  rupture.  This  was  attributed  to  the 
system  being  unable  to  dissipate  the  heat  fast  enough,  resulting 
in  a  high  temperature  and  strain  on  the  foil  due  to  the 
difference  in  the  thermal  expansion  coefficient  of  nickel  and 
tantalum.  A  new  design  was  then  proposed  for  higher  energy 
studies.  It  consisted  of  making  a  hole  in  both  faces  of  the  cell 
and  catching  the  beam  in  a  Faraday  Cup.  Preliminary  work  has 
shown  the  new  design  to  be  successful. 
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CHAPTER  III 

NEUTRON  DETECTOR  EFFICIENCY  MEASUREMENTS 

3.1  The  Relative  Efficiency  Experiment 

Fundamental  to  angular  distribution  studies  is  the 
ability  to  renormalize  the  raw  data  such  that  all  parameters 
giving  rise  to  a  variation  in  count  rate  are  removed,  except 
that  due  to  the  differential  cross  section.  Corrections  for 
variations  in  target  thickness,  charge  collection,  and  run 
time  are  achieved  by  placing  a  secondary  detector  at  a  fixed 
angle  and  normalizing  the  primary  data  to  one  of  the  peaks  in 
the  monitor  spectrum.  The  correction  for  the  variation  in 
detector  efficiency  as  a  function  of  neutron  energy  is 
achieved  by  first  evaluating  the  efficiency  function.  A 
number  of  experimental  techniques  have  been  employed  at  this 
laboratory,  with  varying  degrees  of  success,  in  an  effort  to 
measure  this  function. 

7  7 

The  initial  approach  was  to  use  the  Li  (p,n)Be  reaction 
over  a  proton  energy  range  of  2.5  to  5.0  MeV .  This  gives,  at 
zero  degrees,  a  neutron  energy  range  of  0.786  to  3.331  MeV. 

In  addition  to  the  main  detector  the  neutrons  were  detected  by 
a  rf lat-response"  long  counter  of  the  McKibben  type  (A160) 
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placed  along  side  the  principal  detector.  The  basic  idea  was 

to  unfold  the  efficiency  function  by  normalizing  the  counts  of 

the  main  detector  to  those  of  the  long  counter.  To  make  the 

measurement  meaningful,  however,  one  must  account  for  the  fact 

that  the  long  counter  sees  neutrons  from  both  the  ground  and 

7 

first  excited  states  of  Be  (Be  61).  Whereas,  to  establish  a 
unique  neutron  energy,  only  the  ground  state  peak  of  the  princi¬ 
pal  counter  is  used.  In  addition,  one  must  account  for  the 
known  response  of  the  long  counter  (Al  60) . 

When  a  detailed  evaluation  of  the  efficiency  curve  was 
measured  in  this  way  a  number  of  anomalies,  in  the  form  of 
resonant  humps,  occurred  around  2.5  MeV.  In  addition,  the 
reproducibility  of  points  was  not  very  satisfactory. 

In  an  effort  to  remove  these  anomalies  the  same  technique 
was  applied  to  the  monoenergetic  neutron  sources  of  the 
C12 (He 3,n) O14 ,  C13(a,  n)015,  and  T(p,n)He3  reactions.  However, 
the  cross  sections  for  the  C'*'2  (He3  ,n)  O34  and  C^3  ( a , n) O^3  reac¬ 
tions  proved  to  be  too  small,  leading  to  a  high  background  count 
in  the  long  counter  and  consequently  a  very  poor  reproducibility 
of  points. 

The  tritium  target  was  of  the  zirconium  backed  variety  and, 
due  to  the  energy  loss  of  the  incident  protons  in  penetrating 
the  backing,  gave  rise  to  a  neutron  peak  with  a  long  low-energy 
tail  (Fi  66) .  Specifically,  the  area  of  the  tail  varied  from 
12  to  20  percent  of  the  total  area.  This  then  invalidated  the 


1  I  I 


■ 


■ii  eUm  oT  .*«-*  «  oo  vaoX  «4»  io  »  ort*  toloefb  «x*m  si 


,„i„  ortJ  »o  i  aq  *«•+«  iwotp  sri*  Yi«0  .W  *'  »•  »*» 


■ 


-16- 


idea  of  a  unique  neutron  energy  and  also  made  it  difficult  to 
establish  a  consistant  background  cut-off  criterion. 

The  technique  finally  adopted,  and  that  proved  reasonably 

successful,  consisted  of  an  angular  distribution  study  of  the 

3  3 

T (p ,n) He  reaction  at  3.00  MeV.  and  the  D(d,n)He  reaction  at 

3.52  MeV.  The  efficiency  curve  was  obtained  by  normalizing 

the  experimental  angular  distributions  to  the  published  values 

3 

for  the  respective  reactions  (Go  61),  (Br  60).  The  T(p,n)He 
reaction  covered  a  neutron  energy  range  from  0.728  to  2.62  MeV 

3 

and  the  D(d,n)He  reaction  from  2.31  to  6.79  MeV.  Fig.  3-7 
illustrates  the  relative  efficiency  curve  obtained  in  this  way. 

3 . 2  The  Absolute  Efficiency  Experiment 
( i)  Introduction 

The  basic  inadequacy  in  all  of  the  aforementioned  techniques 
for  measuring  neutron  detection  efficiencies  is  their  dependence 
on  published  data.  To  overcome  this  dependence  and  at  the  same 
time  obtain  absolute  efficiencies  of  high  reliability  a  project 
was  initiated  to  measure  detection  efficiencies  via  the  associ¬ 
ated  particle  technique.  The  object  of  this  technique  is  to 
measure  the  charged  particle  counterparts  of  a  neutron  producing 
reaction  in  coincidence  with  the  neutrons.  From  the  one  to  one 
relationship  that  results  one  can  obtain  absolute  efficiencies 
by  forming  the  ratio  of  the  number  of  neutrons  produced,  as  given 
by  the  charged  particle  counts,  to  the  number  of  neutrons 
detected,  as  given  by  the  neutron  spectrometer. 
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In  choosing  a  reaction  compatible  with  the  demands  of 
this  experiment,  the  following  criteria  were  considered: 

(i)  the  reaction  should  produce  neutrons  within  the 
energy  range  of  interest  (0.5  to  10.0  MeV. ) ; 

(ii)  the  mass  of  the  charged  particle  should  be  small 
so  as  to  share  a  reasonably  high  fraction  of  the 
recoil  energy; 

(iii)  the  reaction  should  be  a  monoenergetic  neutron 
source;  and 

(iv)  to  facilitate  the  escape  of  the  charged  particles 
the  target  material  should  be  in  gaseous  form. 

These  conditions  were  most  successfully  satisfied  by  the 
3 

D(d,n)He  reaction  (Q  =  3.266  MeV.)  and  a  gas  target. 

( ii )  The  Gas  Target  Assembly 

In  addition  to  the  judicious  choice  of  a  reaction, 

the  absolute  efficiency  experiment  required  a  rather  unique 

target  design.  The  principal  requirements  were: 

(i)  facility  for  handling  a  gaseous  target  material 
within  a  vacuum  environment; 

(ii)  containing  the  deuterium  gas  in  such  a  way  and  by 

3 

such  a  material  as  to  allow  the  He  particles  to 
escape ; 

(iii)  provision  for  mounting  and  positioning  a  solid  state 
detector  and  collimation  system;  and 

(iv)  provision  for  aligning  the  target  with  respect  to 
the  incident  beam. 
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In  addition  to  fulfilling  these  requirements,  consideration 
was  given  to  a  scheme  designed  to  maximize  the  active  life 
of  the  target.  This  consisted  of  a  provision  for  rotating  the 
gas  cell,  with  respect  to  the  incident  beam,  and  thus  reducing 
the  weakening  effect  of  the  beam  on  the  very  thin  membrane 
which  served  to  contain  the  gas. 

Fig.  3-1  is  a  relief  drawing  of  the  base  assembly  and  gas 
target  proper.  This  structure  in  turn  fits  up  into  an  existing 
bell  shaped  enclosure  positioned  at  the  end  of  the  Mobley  com¬ 
pression  magnet.  When  so  positioned  an  O-ring  seal  on  the 
upper  face  of  the  bottom  flange  (1)  allows  for  the  subsequent 
evacuation  of  the  target  chamber. 

The  cylindrical  gas  cell  (2)  was  of  the  first  generation 
design.  It  allowed  for  charged  particle  detection  at  all 
angles  except  for  5  degree  cones  at  zero  and  180  degrees.  The 
gas  within  the  cell  was  contained  by  an  ultra-thin  window  of 
Formvar  foil  (3)  which  was  stretched  across  an  annular  gap 
0.200"  high.  The  gap  was  maintained  by  two  tungsten  pegs  (4) 
(0.090"  in  dia.)  positioned  180  degrees  apart.  To  facilitate 
collimation  of  the  incident  beam  and  to  prevent  it  from  hitting 
the  brass  cell,  shields  of  0.010"  gold  (5)  were  placed  around 
the  target  and  slightly  overlapping  the  Formvar  window. 

The  remainder  of  the  target  assembly  consists  of  three 
concentric  shafts,  with  the  inner  two  being  allowed  to  rotate. 


' 


p  a  t  i£  ?•  o:  o 

ao.,.  J  15  o,  i  50  a  M  <£>.  X«  3  »«=  Isoxil.n  I*? 


■ 


INCHES 

1 - 1 - 1 - 1 

0  12  3  4 


DETECTOR® 


FIG.  3- 1:  GAS  TARGET 
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The  function  of  the  center  shaft  (6)  is  to  rotate  the  target 
cell  via  the  variable  speed  motor*  (7)  and  to  provide  a  path 
for  the  entering  gas.  For  smooth  operation  the  shaft  is  moun¬ 
ted  on  two  sets  of  tapered  bearings  (8)  and  (9) . 

The  function  of  the  second  shaft  (10)  is  to  position  and 
monitor  (via  the  vernier  (11))  the  solid  state  detector  (12)  . 

It  will  be  noted  that  two  sets  of  O-rings  (13)  and  (14)  and 
differential  pumping  (15) are  used  in  sealing  these  shafts. 

The  flange  (16) ,  on  the  second  shaft,  and  its  associated 
threaded  mechanism  (17)  allows  one  to  position  simultaneously 
the  target  cell  and  solid  state  detector  in  a  vertical  direction 
and  thus  align  it  with  respect  to  the  incident  beam.  Bearing 
(18)  and  the  teflon  washer  (19)  are  incorporated  to  reduce  the 
operating  friction  of  the  second  shaft. 

As  noted,  the  gas  inlet  (20)  is  provided  with  two  channels: 
one  for  entering  the  gas  (which  is  in  turn  monitored  by  a 
mercury  monometer)  and  another  for  evacuating  the  gas  cell  and 
target  chamber  simultaneously. 

For  the  associated  particle  experiment  a  well  defined  beam 
of  He  particles  is  required.  To  this  end  collimator  (21)  forms 
part  of  the  solid  state  detection  system.  It  is  cone  shaped 
with  two  defining  apertures  (22)  and  (23).  Its  continuous 
stainless  steel  structure  serves  to  differentiate  against  multi- 

Zero-max  motor,  Model  Al-Ml;  Zero-max,  15  Musgrove  St., 
Toronto  13,  Ontario. 
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ply  scattered  particles  and  thus  reduces  the  charged  particle 
background.  In  addition  its  cone  shape  allows  one  to  position 
the  detector  to  within  5  degrees  of  the  beam.  The  two  apertures, 
which  are  easily  replaced,  subtend  an  angle  of  about  4  degrees. 

This  is  a  reasonable  compromise  between  energy  resolution  and 
counting  rate.  The  spacing  and  size  of  the  collimator  aper¬ 
tures  form  the  heart  of  the  geometric  problem  associated  with  the 
efficiency  experiment. 

The  cylindrical  gas  cell  (2)  of  Fig.  3-1  functioned 
reasonably  well,  but  was  subject  to  two  major  weaknesses.  Firstly, 
it  was  no  easy  task  to  produce  a  completely  homogeneous  and  yet 
leak  tight  window  over  the  entire  circumference  of  the  gap.  And 
secondly,  undue  noise  was  created  in  the  solid  state  detector 
when  the  beam  scattered  off  the  tungsten  pegs.  To  overcome  these 
difficulties  the  alternate  design  of  Fig.  3-2  was  used. 

It  consists  of  two  brass  caps  (1)  and  (2)  epoxied*  to  a 
tantalum  cylinder  (3) ,  1"  high  and  1"  in  diameter.  Along  one 
side  of  the  cylinder  is  cut  a  1/4”  wide  grove  over  which  is 
placed  a  Formvar  membrane  (4) .  In  addition  a  1/4"  hole  is  made 
in  the  front  and  back  face  of  the  cylinder.  These  in  turn  are 
covered  by  sheets  of  50  yin.  nickel  foil.  After  passing  through 
a  circular  collimator,  the  incident  beam  enters  and  exits  the 


*Plastic  Steel,  Devcon  Corp. 
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gas  cell  through  these  nickel  windows  and  is  in  turn  stopped 

by  a  Faraday  cup.  The  He  particles,  however,  escape  through 

the  much  thinner  Formvar  window  and  are  in  turn  detected  by 

the  solid  state  system  of  Fig.  3-1. 

Since  the  target  does  not  rotate  the  Far.mvar  foil  is  not 

subjected  to  the  weakening  action  of  tfye  incident  beam.  In 

*■ 

addition  the  area  of  the  Formvar  window  is  substantially  re¬ 
duced.  However,  one  is  restricted  in  the  angular  range  (30 
to  150  degrees)  over  which  the  charged  particles  can  be 
detected . 

(iii)  Construction  of  the  Formvar  Window 

In  preparing  the  Formvar  foil  (Cj-H^C^)  a  procedure  similar 
to  that  outlined  by  Hoogenboom  (Ho  61)  was  adopted.  Fig.  3-3 
is  a  schematic  outline  of  the  apparatus  used. 

First  a  2%  solution  (by  weight)  of  Formvar 15/95*  with 
ethylene  dichloride  is  prepared.  The  solution  is  then  dropped, 
via  a  burette  (20  drops/min.),  onto  a  file  partly  submerged  in 
a  bath  of  distilled  water.  (The  action  of  the  file  being  to 
spread  the  solution  as  it  enters  the  bath.)  Upon  contact  with 
the  water  a  surface  tension  interaction  quickly  spreads  the 
solution  into  a  very  thin  layer,  resulting  in  the  rapid  evapo- 


Obtained  from  Shawinigan  Resins  Corp.,  Springfield,  5, 
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ration  of  the  ethylene  dichloride.  The  resultant  ultra-thin 
(^0.06y)  film  of  about  3  to  4  inches  wide  then  proceeds 
towards  the  rotating  target  cell  as  more  solution  enters  the 
bath.  To  prevent  the  edges  of  the  Formvar  ribbon  from  folding 
back  over  the  window,  as  it  is  being  picked  up  by  the  rotating 
cell,  a  glass  rod  is  introduced  to  guide  the  edges  of  the  foil 
away  from  the  window. 

The  initial  pickup  action  of  the  foil  by  the  rotating 
cell  is  achieved  by  placing  the  surface  of  the  cell  about  1/8" 
below  the  surface  of  the  water.  Once  the  rolling  action 
begins  a  continuous  ribbon  is  easily  achieved  by  the  pulling 
action  of  the  rotating  cell.  However,  the  greatest  difficulty 
is  in  obtaining  a  continuous  foil  of  homogeneous  quality. 
(Homogenity  is  considered  an  important  quality,  not  only  to 
guarantee  a  uniform  energy  loss  but  also  to  obtain  a  maximum 
of  strength  for  a  minimum  of  material.)  The  basis  of  the 
problem  is  that  stress  ripples,  which  become  regions  of  over¬ 
lapping  foil  when  rolled  onto  the  cell,  are  generated  by  the 
localized  pulling  action  of  the  cell.  One  can  reduce  this 
action,  somewhat,  by  gently  manipulating  the  outer  edges  of 
the  foil  with  one's  fingers.  Difficulty  is  sometimes  also 
encountered  in  getting  the  foil  to  spread  initially.  However, 
once  again,  with  the  gentle  manipulation  of  one's  fingers 
one  can  encourage  the  foil  towards  the  rotating  cell. 
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Having  obtained  the  desired  number  of  turns  the  foil 
growth  is  terminated  by  either  removing  the  source  (which 
results  in  an  additional  3  or  4  much  thinner  layers  from  the 
residual  solution  on  the  file)  or  by  placing  the  glass  rod 
in  front  of  the  rotating  cell  and  thereby  breaking  the  foil 
at  the  point  of  contact. 

% 

In  preparing  the  target  cell  for  the  rolling  process 
consideration  had  to  be  given  to  the  sealing  properties  of 
the  Formvar  foil  to  the  gas  cell.  A  number  of  schemes  were 
tried,  but  the  greatest  success  was  attained  by  first  wetting 
the  contact  surface  with  the  2%  Formvar  solution,  followed 
immediately  by  the  rolling  procedure.  In  the  event  of  an 
’’abortive  roll”,  ethylene  dichloride  was  used  to  clean  the 
target  cell. 

In  reality  the  overall  procedure  outlined  above  is  far  from 
perfect,  but  with  practice  a  succcess  rate  of  10-20%  is  obtain¬ 
able.  There  is,  however,  a  slight  correlation  between  the 
success  rate  and  the  procedure  used  in  preparing  the  solution. 
That  is,  clear  solutions  (obtained  by  allowing  the  Formvar 
powder  to  dissolve  slowly  in  the  ethylene  dichloride)  rather 
than  the  milky  solutions  (resulting  from  shaking  the  solution 
to  speed  up  the  dissolving  process)  seem  to  be  a  source  of 


greater  success. 

Using  the  procedure  outlined  above  foils  ranging  in  thick- 
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o  t  q 

ness  from  30  to  300  keV,  with  respect  to  5.15  MeV  Pu  a 
particles,  were  constructed.  It  was  found  that  foils  as  thin 
as  50  keV  could  withstand  a  pressure  gradient  of  2  cm.  of 
Hg.  and  a  beam  current  of  1  yamp  without  undue  strain. 

( iv)  Electronics  of  the  Associated  Particle  Experiment 

Fig.  3-4  is  a  block  diagram  of  the  electronic  configuration 

employed  by  the  associated  particle  experiment.  The  coincidence 

imposed  by  the  system  allows  a  neutron  to  be  detected  if  an  only 

3 

if  at  the  same  time  a  He  particle  is  detected,  whereas,  all 

3 

He  particles  are  detected  independent  of  the  neutrons. 

The  neutrons  are  spectrum  analyzed  by  the  usual  time-of- 
f light  technique,  whereas,  the  charged  particles  are  spectrum 
analyzed  by  a  direct  pulse  height  evaluation  of  the  signals 
coming  from  the  solid  state  detector.  The  neutron  detector  is 
the  same  liquid  scintilator.  Philips  XP  1040  phototube  combi¬ 
nation  discussed  in  Section  2.3.  The  charged  particle  detector 
is  a  Nuclear  Diode  PL-2-35-10  surface  barrier  silicon  detector 
with  a  170y  depletion  depth. 

The  experiment  is  run  with  a  D.C.,  rather  than  pulsed,  beam 
of  incident  deuterons.  The  start  pulse  for  the  time-to-amplitude 
converter  (T.A.C.)  comes  from  the  Ortec  time-pickoff  unit  which 
is  connected  in  series  between  the  solid  state  detector  and 
preamplifier.  The  signal  coming  from  the  detector  passes  through 
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the  primary  winding  of  the  pulse  transformer  and  then  onuo 
the  preamplifier .  The  corresponding  pulse  induced  in  the  secon¬ 
dary,  which  is  subsequently  amplified  and  subjected  mo  a  discrimi¬ 
nation  bias,  serves  as  the  timing  signal.  From  the  preamplifier 
the  original  signal  is  once  again  amplified  and  appropriately 
sniped  ' by  the  Ortec  410  "ultimode  Amplifier)  for  entry  into 
the  pulse  height  analyzer  consisting  of  ADC-A  and  the  or.  line 
SDS  -  920  computer. 

The  coincidence  conditions  imposed  by  the  window  and  dis¬ 
crimination  bias  of  S.C.A  =1  serve  to  single  out  the  He'  peak  and 
to  establish  a  low  energy  discrimination  level.  All  the  pulses 
above  the  discrimination  level,  except  for  those  corresponding 
to  the  He'  peak,  are  analyzed  and  stored  in  the  first  1024 
channels  (0-1023)  of  memory.  Those  corresponding  to  the  He ^ 
peak  are  similarly  analyzed  but  are  routed  irio  the  second  1C 24 
channels  (1024-2047)  of  memory  (for  purposes  of  identifying  the 
window  width) .  Having  once  accumulated  the  respective  daua 
for  a  complete  run  the  original  charged  particle  spectrum,  is 
restored  by  adding  together  the  first  two  1024  channels  of 
memory.  The  resultant  spectrum  is  stored  in  the  third  1C24 
channels  (2048-3071)  of  memory  and  is  used  in  establishing  the 

3 

backa round  count  under  the  He'  peax. 

The  stop  pulse  for  the  T.A.C.,  as  before,  is  taken  i r cm 
the  14th  dynode  of  the  phototube.  The  resultant  time  converted 
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signals  are  subsequently  amplified  by  G.P.A.  #1  and  then  spec¬ 
trum  analysed  and  stored  by  the  T.M.C.  210-1024  pulse  height 
analyzer . 

The  coincidence  conditions  imposed  on  the  neutron  spectrum, 
as  established  by  Triple  Coincidence  Unit-B,  are: 

(i)  that  a  signal  fall  within  the  window  setting  of 

S.C.A.  #1  (corresponding  to  the  detection  of  a 
3 

He  particle) ,  and 

(ii)  that  the  photo-election  signal,  taken  from  the  11th 

dynode  of  the  phototube,  be  above  the  discrimination 
level  established  by  S.C.A.  #2. 

The  resultant  neutron  spectrum  is  then  a  single  peak  corres¬ 
ponding  to  those  neutrons  produced  in  coincidence  with  a  kine- 

3 

matically  associated  He  particle. 

At  the  end  of  data  collection  the  neutron  spectrum  stored 
within  the  T.M.C.  210-1024  unit  is  transferred  into  the  fourth 
1024  channel  memory  block  of  the  SDS-920  computer.  Fig.  3-5 
shows  the  essential  features  of  the  spectral  contents  of 
memory.  From  the  computer  memory  the  data  is  transferred  to 
magnetic  tape  for  further  analysis. 

It  is  important  to  the  correct  evaluation  of  the  detector 
efficiency  to  have  the  photo-electron  threshold  level  of 
S.C.A.  #2  equal  to  that  used  by  the  spectrometer  configuration 
of  Section  2.3.  This  criterion  is  fulfilled  by  maintaining  a 


fixed  setting  for  a  parameter  whose  measurement  is  sensitive  to 
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the  discrimination  level  of  S.C.A.  #2.  The  parameter  so 

2  2 

chosen  is  known  as  the  "Na  cut  off  ratio"  and  is  determined 
in  the  following  way. 

The  inputs  of  ADC-A  and  the  output  of  S.C.A.  #2  are  removed 

and  replaced  by  a  parallel  connection  from  the  output  of  G.P.A. 

#2  to  the  analog  input  of  ADC-A  and  a  connection  from  the  out- 

22 

put  of  S.C.A.  #2  to  the  coincidence  input  of  ADC-A.  A  Na  y 

source  (0.511  MeV  from  positron  annihilation  and  1.274  MeV  from 
22 

Ne  decay)  is  then  placed  in  front  of  the  neutron  detector  and 

22 

the  associated  pulse  height  is  collected.  The  ''Na  cut  off 
ratio'* * 3 ,  R^a22  >  is  given  by 


RNa22 


c  -  b 
b  -  a' 


(3.2-1) 


where  a,  b,  and  c  are  the  channel  positions  outlined  in  Fig. 

3-6. 

2  2  . 

A  "Na  cut  off  ratio"  of  3.249  was  used  throughout,  giving 

a  neutron-detection  threshold  energy  of  0.80  MeV. 


3 . 3  Presentation  of  Results 

Having  established  the  desired  coincidence  conditions  a 

simple  minded  approach  to  the  efficiency  measurement  consists 

of  positioning  the  neutron  detector  such  that  it  subtends  all 

of  the  neutrons  produced  in  coincidence  with  the  associated 

3 

cone  of  He  particles. 


The  absolute  efficiency  is  given  by  the 
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number  of  neutrons  detected,  divided  by  the  number  of  He  parti¬ 
cles  detected  (which  corresponds  to  the  number  of  neutrons  inci¬ 
dent  on  the  detector).  Unfortunately,  such  an  approach  is 
subject  to  an  over  establishment  of  the  detector  efficiency  for 
the  following  reasons. 

(i)  To  guarantee  that  all  the  coincident  neutrons  are 

incident  on  the  detector  the  perimeter  of  the  coinci¬ 
dent  neutron  flux  must  be  enclosed  well  within  the 
perimeter  of  the  detector  face.  As  a  consequence 
one  neglects  the  reduced  detection  efficiency  associ¬ 
ated  with  neutron  out  scattering  at  the  edge  of  the 
scintillator . 

(ii)  By  employing  a  target  of  finite  volume  and  a  colli¬ 
mator  of  double  aperture  the  coincident  neutron  flux 
is  not  of  constant  value.  Rather,  it  is  in  the  form 
of  a  flat-topped,  2-dimensional  distribution  of  ellip¬ 
tical  cross  section  (see  Figs.  A- 3  and  A-4) .  Such  a 
distribution  serves  to  complement  the  edge  effect 
error  of  (i)  due  to  a  further  concentration  of  neutrons 
at  the  center  of  the  scintillator. 

(iii)  It  has  been  found  (Me  67)  that  the  response  of  the  XP 
1040  photo  tube  is  not  of  constant  value  with  respect 
to  radial  distance,  but  is  rather  characterized  by 
peaks  and  valleys  with  a  sharp  drop  off  near  the 
edges.  Since  we  are  dealing  with  a  concentration  of 
neutrons  at  the  center  of  the  scintillator,  the  non¬ 
uniformity  of  the  phototube  response  would  tend  to 
enhance  the  detector  efficiency. 


To  eliminate  these  errors  it  becomes  necessary  to  lengthen 
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the  neutron  flight  path  so  that  the  detector  face  is  enclosed 
within  the  flat-topped  portion  of  the  neutron  distribution. 

In  addition  to  having  to  contend  with  a  reduced  counting  rate, 
it  then  becomes  necessary  to  evaluate  a  correction  factor,  R, 
giving  the  fraction  of  the  distribution  subtended  by  the  detector. 
The  absolute  efficiency  is  then  given  by 


_  1_  No.  of  neutrons  detected 


(3.3  -  1) 


R  No.  of  He  particles  detected 
For  the  results  discussed  here  the  parameter,  R,  was  calcu¬ 
lated  via  a  computer  program;  written  by  J.  F.  Easton  and  dis¬ 
cussed  under  Appendix  A. 

Fig.  3-7  gives  the  results  of  the  relative  efficiency 

.  + 
measurement  normalized  to  an  absolute  value  of  12.8%  (-  7%)  at 

2.82  meV.  The  target  geometry,  flight  path,  and  R  value  rela¬ 
tive  to  this  measurement  are  given  on  the  data  sheet,  Fig. 

A-2.  The  corresponding  neutron  and  charged  particle  spectra 
are  those  of  Fig.  3-5. 

In  addition  to  the  experimental  data,  Fig.  3-7  gives  the 
results  of  computer  generated,  theoretical  calculations  (Gr  67). 
The  model  used  accounts  for  multiple  neutron  scattering  within 
the  scintillator  plus  the  influence  of  neutron-carbon  scattering. 
Due  to  the  rather  questionable  agreement  the  experimental  curve 
was  used  throughout  in  the  analysis  of  the  angular  distribution 
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To  check  the  influence  of  points  (i) ,  (ii) ,  and  (iii) 
the  associated  particle  experiment  was  repeated  with  the 
neutron  flux  enclosed  within  the  detector  face.  The  results 
of  this  measurement  gave  an  efficiency  value  of  18.6%  at  2.82 
MeV  as  compared  to  the  accepted  value  of  12.8%. 

3 . 4  Conclusions 

The  results  of  this  work  indicate  that  the  ideas  of  an 
associated  particle  experiment  are  successful  in  obtaining  abso¬ 
lute  neutron  detection  efficiencies  to  within  an  accuracy  of  7%. 
However,  the  technique  is  far  from  perfected  and  the  following 
suggestions  are  given  as  possible  improvements. 

(i)  With  reference  to  Fig.  3-5  we  see  that  the  charged 

particle  background  is  excessively  high,  leading  to 

an  unnecessarily  high  error  in  the  determination  of 
3 

the  He  count.  The  primary  source  of  this  background 
is  believed  to  come  from  multiple  elastic  scattering 
of  the  incident  deuteron  beam  by  the  target  and 
collimation  system.  Thus  if  one  were  to  use  a  solid 
state  detector  with  a  depletion  layer  sufficiently 

3 

thin  to  stop  only  the  He  particles  and  not  the 
deuterons,  the  background  should  be  reduced  signifi¬ 
cantly  . 

(ii)  The  most  glaring  inadequacy  of  the  results  of  Fig.  3-7 
is  that  only  one  absolute  efficiency  point  is  given. 
This  is  a  result  of  the  limited  running  time  allowed 
for  the  experiment  and  the  restrictions  imposed  on 

3 

angle  and  energy  for  which  the  He  peak  was  reasonably 
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isolated.  The  basis  of  this  problem  was  a  small 

concentration  of  oxygen  in  the  deuterium  gas .  The 

resultant  deuteron-oxygen  elastic  scattering  produced 

3 

a  peak  near  the  He  and  tritium  peaks  and  equal  in 

magnitude.  It  was  only  through  a  judicious  choice  of 

3 

angles  that  the  He  peak  could  be  separated  from  the 
scattered  D-D  peak  on  the  one  hand  and  the  scattered 
D-Oxygen  peak  on  the  other.  For  the  angle  chosen 

(31.2°)  the  two  scattered  peaks  merged,  leaving  the 

3  .  . 

He  and  tritium  peaks  reasonably  well  isolated.  To 
overcome  this  problem  a  purer  supply  of  deuterium  is 
suggested . 
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CHAPTER  IV 
2  0  2  1 

THE  Ne  (d,n)Na  REACTION 

4 . 1  Introduction 

The  study  of  the  Ne2 0 (d ,n) Na2 1  reaction  at  Ed  =  3.924  MeV 
was  undertaken  with  the  following  objectives  in  mind; 

(i)  to  initiate  a  search  for  new  low  lying 
levels  in  the  Na^  1  nucleus; 

(ii)  to  obtain  spectroscopic  information 
about  the  Na2 1  nucleus  from  a 
D.  W.  B.  A.  analysis  of  absolute 
cross  section  data;  and 

(iii)  to  corroborate  the  experimental 
results  within  the  framework 
of  the  collective  model. 

The  established  structure  of  Na2 1 ,  and  its  minor  nucleus 
Ne2 1 ,  is  illustrated  in  Fig.  4-1.  The  3/2+  assignment,  for 
the  Na2 1  ground  state,  has  been  established  by  the  superallowed 
character  of  its  £+  decay  to  the  ground  state  of  Ne20(En  62) 
and  also  by  an  atomic-beam,  magnetic  resonance  experiment 
(Am  65).  The  early  results  of  Ne " 0  (d ,n) Na2 1  stripping  experi¬ 
ments  restricted  the  J1  value  of  the  Na2 1  0.335  level  to 
either  3/2+  or  5/2+.  C.  Van  der  Leun  and  W.  L.  Monton  (Va  64) 
singled  out  the  5/2+  value,  via  a  Ne^p ,y ) Na2 1  experiment. 
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The  7/2  value  for  the  Na2 1  1.751  level  has  been  assigned  on 
the  basis  of  mirror-pair  arguments  and  a  Ne 2  0  (p , y)  Na2  1  corre¬ 
lation  experiment  (A1  65).  The  Na2 1  2.433  level  has  been 
readily  identified  from  Ne2  0 (d,n)Na2 1  stripping  experiments 
as  a  &p=0  transfer ,  and  thus  the  l/2+  assignment. 

The  question  of  additional  low  lying  structure  for  the 
Na2 1  nucleus  has  been  suggested  by  the  following  work. 

(i)  J.  B.  Marion  (Ma  56),  in  a  Ne 2 0 (d , n) Na 2 1  threshold 
experiment,  claims  to  have  measured  a  level  at  1.46  -  0.04  MeV. 
He  does  not  claim,  however,  to  have  measured  the  1.75  MeV. 
level  established  by  subsequent  Ne2 0 (d ,n) Na2 1  work  (Aj  61  and 
Be  61) .  R.  E.  Benenson  (Be  61) ,  in  addition  to  studying  the 

Ne 2 0 (d ,n) Na2 1  reaction,  also  studied  the  Ne 2 0 (p , y ) Na2 1  reaction 
and  claims  to  have  measured  a  1.42  -  0.04  MeV.  y-ray.  As  its 
possible  source  he  suggests  either  a  ground  state  decay  from 
the  questionable  1.46  level  or  an  element  of  a  triple  cascade 
from  the  1.75  MeV.  level  to  the  0.33  MeV.  level  after  having 
been  fed  from  the  3.56  MeV.  level.  Although  the  evidence  is 
not  overwhelming,  he  suggests  the  latter  as  the  more  acceptable 
scheme.  This  same  conclusion  is  reached  by  Arnell  and  Sterner 
(Ar  63)  after  repeating  the  Ne 2 0 (p , y) Na 2 1  experiment.  The 
arguments  are  not  conclusive,  however,  and  could  be  disregarded 
if  subsequent  evidence  for  a  Na2 1  1.46  MeV.  level  was  found. 

(ii)  W.  R.  Gibbs  (Gi  65)  in  a  study  of  the  Ne 2 0 (d , ny) Na2 1 
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reaction  suggests  the  slight  possibility  that  either  the 
ground  or  first  excited  state  of  Na2 1  might  be  a  double  state. 

(iii)  F.  Everling  (Ev  63)  in  a  study  of  the  total  nuclear 
binding  of  excited  states  for  the  2s  -  Id  subshells  suggests 
the  possibility  of  a  l/2+  level  at  2.18  MeV.  in  Na2 1 . 

(iv)  Recently  it  has  been  shown  (Po  67  and  Ho  67)  that 
the  Ne21  2.80  MeV.  level  is  in  reality  a  2.797  and  2.790  MeV. 
doublet.  This  evidence  then  supports  possible  doublet  structure 
for  the  2.433  MeV.  level  of  Na2 1 ,  since  it  is  the  mirror-pair 
analog  of  the  Ne2 0  2.80  level. 

Since  the  energy  resolution  capability  of  the  neutron 
spectrometer  used  at  this  laboratory  is  significantly  superior 
to  that  obtained  in  any  previous  Ne20  (d,n)Na21  work,  it  was 
felt  that  further  study  of  this  reaction  might  unearth  possible 
doublet  structure.  The  choice  of  reaction  energy  (E^=3.924  MeV.) 
capitalized  on  the  increased  resolution  associated  with  lower 
energies  and  yet  allowed  for  the  study  of  all  the  bound  states 
of  Na2 1 . 

In  addition  to  a  search  for  new  low  lying  structure,  the 
present  experiment  served  to  provide  additional  absolute  cross 
section  data  to  add  to  the  limited  amount  reported  to  date 
(Aj  61)  and  (Pe  65) .  And  with  the  aid  of  a  Distorted  Wave  Born 
Approximation  (D.  W.  B.  A.)  program,  written  by  B .  E.  F. 
Macefield  (Ma  64) ,  spin  and  parity  assignments  were  checked  and 
absolute  spectroscopic  factors  obtained. 
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To  complete  the  treatment  of  the  Ne2 0  (d ,n) Na2 1  experiment, 
a  discussion  of  the  structure  of  Na21  is  made  within  the  frame¬ 
work  of  the  collective  model.  In  this  context  the  absolute 
spectroscopic  factors  are  interpreted  in  terms  of  the  expan¬ 
sion  coefficients  of  the  Nilsson  model  wave  functions. 

4 . 2  The  Experiment 

The  Ne2 0  (d ,n) Na2 1  experiment,  at  a  reaction  energy  of 
3.924  -  0.01  MeV. ,  was  performed  using  the  time-of-f light  system 
and  gas  target  cell  described  in  Chapter  II.  The  average 
deuteron  beam  current  ranged  from  1.25  to  2.00  yamps.  at  an 
energy  of  4.004  -  0.006  MeV.  (An  energy  loss  of  80  KeV.  could 
then  be  attributed  to  the  50  yin  nickel  foil  and  neon  gas.) 

The  Ne20  target  consisted  of  the  gas  cell  (Fig.  2-3) 
filled  to  a  pressure  of  about  43  cm.  of  Hg.  (at  27°C)  with 
99.9  molar  %  Ne2  .  This  corresponded  to  a  target  thickness 
of  about  290  ygm/cm2 . 

The  neutrons  were  timed  over  a  flight  path  of  6.05  meters 
at  19  different  angles  ranging  from  zero  to  140  degrees.  For 
each  angle  a  charge  of  1500  ycoul.  was  collected.  A  timing 
resolution  of  0.96  ns.  for  the  4.01  MeV.  zero-degree,  ground- 
state  neutrons  was  obtained. 

The  monitor  detector  was  positioned  at  -20°  with  respect 
to  the  incident  beam  and  at  a  distance  of  3  meters  from  the 

k 

Obtained  from  Monsanto  Research  Corp. ,  Mound  Laboratory, 
Miamisburg,  Ohio,  U.  S.  A. 
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target.  Due  to  its  strength  and  isolation,  the  Na21  2.433  MeV. 
monitor  peak  area  was  used  as  the  angular  distribution  normali¬ 
zation  factor. 

A  calibration  of  the  kicksorter  was  achieved  by  running  the 
Be9(d,n)B10  reaction,  at  3.00  MeV,  immediately  following  the 
Ne 2 0 (d, n) Na 2 1  experiment.  The  flight  times  for  the  respective 
levels  of  B10  were  calculated  by  a  relativistic  kinematic  program 
(Gr  67a)  on  the  basis  of  the  measured  neutron  flight  path  (6.05  M.) 
and  the  known  level  structure  of  B10  (Aj  66).  The  corresponding 
peak  positions  were  then  least  squares  fitted,  by  a  general  least 
squares  fitting  program  (Da  66) ,  to  a  quadractic  equation  relating 
channel  position  to  flight  time. 

The  Ne2 0 (d,n) Na2 1  reaction  energy  of  3.924-0.01  MeV.  was 
determined  by  first  using  the  calibration  equation  to  obtain  the 
flight  time  corresponding  to  the  ground  state  peak  of  the  run  taken 
at  20  degrees.  A  variation  in  the  reaction  energy  required  by  the 
kinematic  program  was  then  made  until  the  kinematically  generated 
flight  time  was  equal  to  that  determined  by  the  calibration 
equation.  The  ground  state  was  used  in  this  determination  since, 
of  the  levels  excited,  its  Q  value  is  known  to  the  highest  degree 
of  accuracy  (0.208-0.008  MeV.  (La  66)). 

Fig.  4-2  is  the  neutron  spectrum  corresponding  to  a  detec¬ 
tion  angle  of  20  degrees,  and  is  characteristic  of  the  spectra 
obtained.  The  contaminate  levels  due  to  oxygen  and  carbon  are 
believed  to  come  from  either  the  nickel  foil  or  tantalum  backing 
used  in  the  construction  of  the  gas  cell.  It  should  be  pointed 
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out  that  a  parallel  time-of-f light  spectrum  was  formed  in  the 
first  1024  channels  of  memory.  It  consisted  of  those  signals 
incorrectly  identified  by  the  n  -  y  discriminator  as  y-rays. 

To  within  the  resolution  indicated  by  Fig.  4-2,  there  is 
no  evidence  to  support  the  existence  of  doublet  structure  in 
either  the  ground,  first,  or  third  excited  states.  In  addition 
there  is  no  evidence  to  support  the  questionable  1.42  or  2.18 
MeV.  levels.  These  conclusions  are  borne  out  by  the  spectra 
obtained  at  all  angles. 

A  peak  fitting  program  developed  by  J.  W.  Tepel  (Te  66 
and  Gr  66)  was  used  in  determining  the  peak  areas  and  peak 
positions  of  the  time-of-f light  spectra.  The  program  is  so 
designed  such  that  the  shape  of  the  analytic  function  used 

in  the  fitting  process  is  determined  by  a  characteristic  refer¬ 
ence  peak.  For  this  the  Na  *  1  2.4  33  peak  of  the  5  degree  run  was 
used.  Consistent  and  objective  fits  were  obtained  by  using 
this  same  reference  peak  throughout. 

Peak  positions  obtained  in  this  way  were  then  converted 
to  flight  times  via  the  calibration  equation  and  subsequently 
used  by  a  least  squares  mass  and  Q-value  program  (Da  66)  to 
obtain  the  mass  and  excitation  energy  of  the  residual  nucleus. 
The  Q-values  and  excitation  energies,  obtained  in  this  way, 
are  listed  in  Table  4-1  along  with  previously  determined  values. 
(It  should  be  noted  that  the  excitation  energies  listed  are 
based  on  the  known  mass  value  of  the  residual  nucleus.) 
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TABLE  4-1 


Q  Values  and  Excitation  Energies 


Nuclear 

THIS  EXPERIMENT 

PREVIOUS  WORK 

Level 

Q  Value 
(MeV) 

Excitation 

(MeV) 

Q  Value 
(MeV) 

Excitation 

(MeV) 

Reference 

gnd 

0. 

208-0.008* 

0.230-0.03 

Aj 

61 

0.220-0.03 

Be 

61 

1st 

-0. 

127-0.01 

0.335-0.01 

-0.140-0.04 

0.370-0.04 

Aj 

61 

-0.080-0.03 

0.330-0.03 

Be 

61 

0.338-0.003 

Va 

64 

0.342-0.008 

Gi 

65 

2nd 

-1. 

543-0.015 

1.751-0.015 

-1.460-0.04 

1.730-0.04 

Aj 

61 

-1.550-0.05 

1.770-0.04 

Be 

61 

3rd 

-2. 

225-0.01 

2.433-0.01 

-2 .201-2 .007 

2.431-0.031 

Aj 

61 

-2 . 200-0 . 04 

2.420-0.04 

Be 

61 

*This  value  was  not  determined  in  this  experiment  but  is  that 
given  by  (La  66) .  All  excitation  energies  listed  are  based 
on  this  ground  state  Q  value. 
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Absolute  cross  section,  angular  distribution  data  was 
obtained  for  the  Ne2 0 (d,n) Na2 1  reaction  by  subjecting  the  peak 
fitted  results  to  the  following  corrections  : 

(i)  the  neutrons  incorrectly  identified 
as  y-rays  were  added  to  their 
neutron  spectrum  counterparts; 

(ii)  the  resultant  sum  was  divided  by 

the  area  of  the  Na2 1  2.433  monitor 
peak.  (This  corrected  for  any 
variation  in  target  thickness, 
charge  collection,  or  run  time.); 

(iii)  a  correction  due  to  the  difference 
in  live  times  between  the  main 
and  monitor  systems  was  made ; 

(iv)  a  correction  due  to  a  variation 

in  detection  efficiency  was  made ; 

(v)  and  finally,  the  resultant  relative 
angular  distributions  were  normal¬ 
ized  to  an  absolute  determination 
made  at  20  degrees. 

Tables  4-2  to  4-5  list  the  results  of  this  work. 

Most  reactions  designated  as  "direct'’  have  a  small  compon¬ 
ent  of  compound  reaction  formation  present.  This  compound  reac¬ 
tion  component  is  energy  dependent  as  identified  by  strong  reso: - 
ant  structure,  in  the  absorption  cross  section,  around  energies 
corresponding  to  the  population  of  excited  states  in  the  compound 
nucleus.  As  the  absorption  cross  section  goes  up  the  cross 
section  in  the  exit  channels  also  goes  up.  Thus  by 
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TABLE  4-2 


Ne20 (d,n)Na^  Ground  State 


0(Lab)  Deg 

dcr/d!7  rob 
(Lab)  St 

0(c.m. )  Deg 

?a/d.f‘  mb 
(c.m) 

da/d^ 

%  error 

0.27 

1.43 

0.29 

1.25 

10.0 

5.03 

- 

5.37 

— 

- 

1 0  •.  2  8 

1.38 

10.97 

1.21 

8.2 

15.11 

1.27 

16.24 

1.12 

9.4 

20.73 

1.15 

22.11 

1.02 

8.5 

25.82 

1.07 

27.51 

0.951 

8.2 

30.20 

0.919 

32.16 

0.820 

9.5 

36.02 

0.895 

38.31 

0.805 

9.0 

41.20 

0.833 

43.76 

0.755 

9.5 

45.40 

0.767 

48.17 

0.700 

10.1 

51.27 

0.635 

54.30 

0.585 

10.4 

60.75 

0.545 

64.14 

0.511 

12.8 

71.45 

0.521 

75.14 

0.500 

13.1 

75.65 

0.524 

79.42 

0.508 

12.6 

90.55 

0.564 

94.44 

0.566 

11.0 

100.43 

0.628 

104.26 

0.645 

9.9 

115.45 

0.637 

118.96 

0.677 

9.8 

120.00 

0.699 

123.37 

0.75C 

10.0 

139.53 

0.579 

142.05 

0.644 

- - 

10.4 

~0  £*  — 
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j.  (ra.a) 
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rs.o 

«pa 

Y£.2 
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XS.X 
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u.ss 

rr.os: 
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TABLE  4-3 


Ne2  0  (d  ,n)  Na2  1_ 0.33  5  MeV  Level 


t 

0 (Lab)  Deg 

do  /dQ 
,T  i  \  mb 

(Lab)  St 

0(c.m.)  Deg 

mb 

(c.m.)  ^ 

do/dQ 
%  Error 

0.27 

1.39 

0.29 

1.21 

11.9 

5.03 

1.44 

5.39 

1.26 

10.3 

10.28 

1.75 

11.01 

1.53 

10.0 

15.22 

2.37 

16.29 

2 . 03 

9.9 

20.73 

3 . 60 

22.17 

3.17 

9.4 

25.82 

4.82 

27.59 

4.26 

8.7 

30.20 

5.73 

32.25 

5.09 

8.5 

36.02 

6.36 

3  8r.  4  2 

5.69 

7.9 

41.20 

6.69 

43.88 

6.03 

7.7 

45.4  0 

6.50 

48.30 

5.90 

8.6 

51.27 

5.99 

54.45 

5.50 

9.1 

60.75 

4.72 

64.30 

4.42 

9.0 

71.45 

3.19 

75.31 

3.06 

9.4 

75.65 

2.43 

79.60 

2.35 

8.3 

90.55 

1.49 

94.62 

1.50 

9.0 

100.43 

1.20 

104.44 

1.23 

9.0 

115.45 

1.11 

119.13 

1.18 

8.9 

120.00 

1.23 

123 . 53 

1.326 

9.1 

139.53 

1.14 

142.17 

1.27 

9.6 
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TABLE  4-4 


Ne20  (d  ,n)  Na2  1  1.751  MeV.  Level 


0(Lab)  Deg 

da/ b 
(Lab,  || 

e(c.m.)  Deg 

?£/dn.  mb 

da/d  ft 
%  Error 

0.27 

0.439 

0.30 

0.368 

12.8 

5.03 

0.395 

5.50 

0.331 

14 . 8 

10.28 

0.454 

11.23 

0.381 

13.8 

15.22 

0.478 

16.61 

0.403 

11.8 

20.73 

0.400 

22.61 

0.339 

12.4 

25.82 

0.496 

28.13 

0.423 

10.3 

30.20 

0.570 

32.87 

0.489 

11.9 

36.02 

0.552 

39.14 

0.478 

17.0 

41.20 

0.557 

44.70 

0.488 

15.0 

4  5.40 

0.512 

49.18 

0.452 

13.5 

51.27 

0.565 

55.41 

0.506 

12.8 

60.75 

0.392 

65.38 

0.360 

15.1 

71.45 

0.460 

76.48 

0.436 

17.7 

| 

75.65 

0.511 

80.79 

0.490 

16.5 

90.55 

0.560 

95.86 

0.563 

12.8 

100.43 

0.456 

105.65 

0.474 

13.5 

115.45 

0.345 

120.24 

0.376 

14.3 

120.00 

0.366 

124.60 

0.404 

14.4 

139.53 

0.253 

142.97 

0.293 

22.4 
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TABLE  4-5 


Nez  0  (d  ,n)  Na2  x _ 2.433  MeV.  Level 


0 (Lab)  Deg 

da/d^  mb 

(Lab)  St 

(c.m.)  Deg 

dcr/diw!  mb 

(c.m.)  St 

da/d^ 

%  Error 

0.27 

61.85 

0.30 

49.87 

7.1 

5.03 

57.32 

5.60 

46.26 

7.2 

10.28 

47.52 

11.44 

38.45 

7.3 

15.22 

39.04 

16.93 

31.72 

7.2 

20.73 

25.29 

23.04 

20.68 

7.2 

25.82 

16.08 

28.66 

13.25 

6 . 6 

30.20 

11.29 

33.48 

9.37 

7.3 

36.02 

6.78 

39.85 

5.69 

8 . 5 

41.20 

5.88 

45.49 

5.01 

11.2 

45.40 

5.92 

50.04 

5.09 

11.8 

51.27 

7.10 

56.36 

6.21 

9.6 

60.75 

8.54 

66.44 

7.70 

9.7 

71.45 

8.26 

77.64 

7.73 

8 . 2 

75.65 

7.66 

81.97 

7.28 

7 . 6 

90.55 

5.36 

97.07 

5.40 

9.4 

100.43 

4.11 

106.85 

4.31 

10.5 

115.45 

2.34 

121.38 

2.60 

8.4 

120.00 

2.43 

125.65 

2.75 

8 . 6 

139.53 

1.65 

143.76 

1.98 

9.5 

32.0 
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measuring  the  Ne 2 0  (d ,n) Naz 1  cross  section  with  respect  to  a 
variation  in  reaction  energy  one  can  identify  the  presence 
of  strong  compound  reaction  formation.  Since  in  an  effective 
D.  W.  B.  A.  analysis  one  wishes  to  minimize  the  presence  of 
compound  reaction  formation  one  searches  for  a  minimum  in  the 
Ne 2 0 (d , n) Na 2 1  cross  section  with  respect  to  energy.  To  this 
end  a  single  angle  excitation  experiment  was  performed.  The 
energy  range  =  3.856  to  4.016  MeV.  was  covered  at  a  detec¬ 
tion  angle  of  20.73  degrees.  Fig.  4-3  gives  the  resultant 
curves  for  the  ground,  first,  and  third  excited  states. 

Although  total  cross  section  data  is  required  before  a  firm 
conclusion  can  be  drawn,  the  results  of  Fig.  4-3  tend  to 
support  the  assumption  of  a  constant  background  contribution. 

4 . 3  The  D.  W.  B.  A.  Analysis  of  the  Ne 2 0  (d ,n) Na 2 1  Reaction 

(i)  D.  W.  B.  A.  Theory 

The  distorted-wave  Born-approximation  (D.  W.  B.  A.)  theory 
is  a  description  of  a  class  of  reactions  known  as  "direct". 

The  interaction  responsible  for  the  reaction  is  assumed  to 
occur  just  once,  with  the  transition  taken  between  states  of 
elastic  scattering.  The  theory  then  corresponds  to  a  one  step 
process . 

The  formal  expression  for  the  transition  complitude  corres¬ 
ponding  to  the  generalized  reaction,  A(a,b)B,  as  developed  by 
R.  H.  Bassel  (Ba  64)  is: 


6  39'Jqft' n  :‘^iv  n  )i^oa«  eec  io  £  1  1  e  9  UB:  71 
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TDW  -  J/di/d?b  Xb(’)(5b'  ?b>  “  Bb|v|Aa>x<  +  )  (ka,  ra)  (4.3-1) 


Here  r&  is  the  displacement  of  'a1  from  'A'  and  rfe  the  displace¬ 
ment  of  '  b'  from  '  B'.  J  is  the  Jacobian  of  the  transformation 

to  these  relative  coordinates.  The  functions  y  ^ ^  and  y/  ^ 

b 

are  the  elastic  scattering  states,  where  (+)  and  (-)  denote 
the  usual  outgoing  and  incoming  boundary  conditions.  V  is  the 
interaction  potential  and  | Aa >  and  | Bb >  are  the  internal  states 
of  the  colliding  pairs  before  and  after  the  reaction. 

In  the  optical-model  approximation  the  scattering  state 
X^+)  and  ^  are  solutions  of  the  Schrodinger  equation: 

[v2+k2-  (U(r)  +  u(r)  )  ]  X  (K,?)  =  0  (4.3-2) 


where  U(r)  is  the  optical-model  potential,  (r)  the  coulomb 
potential,  and  y  the  reduced  mass  of  the  pair. 

In  its  most  generalized  form  the  optical-model  potential 
is  given  by: 


U (r ) =  - [Vf (r)  +  i  Wg(r)  +  4iW  g '  (r)  +  X  (S -L)  df  ]  (4.3-3) 

D  dr 

where  f(r),  g(r),  and  g'(r)  are  the  Wood-Saxon  form  factors 
given  by: 
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g'  (r) 


xw 

(1+e  *) 


2 


(4.3-6) 


The  first  tern  of  eqn.  (4.3-3)  corresponds  to  volume  potential 
scattering,  the  second  to  volume  potential  absorption,  the 
third  to  surface  absorption  and  the  fourth  to  spin-orbit  scat¬ 
tering  . 

2  0  2  1 

In  the  treatment  of  the  Ne  (d,n)Na  data  the  terms 
W  and  X  were  set  to  zero,  consistent  with  the  ideas  of  a  purely 
surface  absorption  interaction  and  negligible  spin-orbit  effect. 

By  setting  X  (which  in  general  is  complex)  equal  to  zero 
one  can  neglect  the  spin  structure  of  the  particles.  The 
resulting  wave  function,  x_  V  ^  and  Xu  ^  ^  '  are  then  reduced  from 
matrices  to  scalar  functions  and  are  related  in  the  following 
way  by  the  condition  of  time  reversibility. 


xl_)*<5,5)  =  x<  +  )  (-5,5). 


(4.3-7) 


In  performing  the  integration  of  eqn.  (4.3-1)  the  scatter¬ 
ing  state  wave  functions  are  first  expanded  into  a  series  of 
partial  waves,  i.e. : 


X 


(  +  ) 


i!  I  iL 

k  L,M 


XL(k,r) 


VM<£, 


Yl  <?) 


(4.3-8) 


The  resultant  radial  wave  function,  x_  (k,r) ,  is  given  by: 


il  +  k"  -  L(L~t'1)-  -  — (V(r)  +  iW  (r )  +  U  (r) 
dr2  r2  **  Dc 


XT  (k,r) =0  (4 . 3-9) 


and  takes  on  the  asymtotic  form: 


-47 


XL(k,r)  =  j  [HL*(kr)  -  r^H^kr)]  e1(?L 


(4.3-10) 


where  H7  =  (G  +  i  F  ) ,  the  Coulomb  functions  irregular  and 

regular  at  the  origin;  a  is  the  Coulomb  phase  shift;  and  n  is 

L  L 

the  reflection  coefficient  for  the  L1^  partial  wave. 

The  remaining  factor  in  eqn.  (4.3-1)  is  the  interaction 

matrix  element,  <Bb|v|Aa>.  This  factor  is  a  function  of  r  and 

a 

and  plays  the  role  of  an  effective  interaction  for  the  transi¬ 


tion  between  the  elastic  scattering  states,  xa  an<^  Xb  •  It 
contains  all  the  information  on  nuclear  structure,  angular  momen¬ 
tum  selection  rules  and  even  the  the  type  of  reaction  being  con¬ 
sidered  (i.e.  stripping,  pickup,  or  inelastic  scattering). 

A  number  of  general  statements  can  be  made  about  the  matrix 
element  by  expanding  it  into  terms  where  a  definite  angular 
momentum,  ] ,  is  transferred  to  nucleus  A  to  form  the  residual 
nucleus  B.  That  is 


5  =  J. 


J 


B 


(4 »  3-11) 


where  j  is  made  up  of  an  orbital  angular  momentum  transfer  l , 

— ^  c 

and  a  spin  part  s,  given  by; 


s  =  s  -  s, 
a  b 


(4.3-12) 
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l  =  j  -  s 


(4.3-13) 
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The  product  J<Bb|v|Aa  is  then  expressed  as: 
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J<sn 

m,  , 
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j  M 
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’G£sj,m(5b'  V  bB'  aA) 


B  A 


(4.3-14) 


where  the  Clebsch-Gordon  coefficients  take  care  of  the  angular 

momentum  coupling ,  and  where  m=M  -M  +m  -  m,  .  The  notation, 

^  ^  B  A  a  b 

bB  and  aA,  denotes  the  possible  dependence  of  G  on  other  nuclear 

-  £ 

and  projectile  quantum  numbers  and  the  factor  i  is  included  to 
ensure  convenient  time-reversal  properties.  Formally  Q  is 

36  S  "J  f  m 

defined  by  the  inversion  of  eqn.  (4.3-14)  (Sa  64). 

In  general  the  function  G„  .  can  be  factored  into  an 

£s;j  ,m 

amplitude  and  a  form  factor: 


G£sj ,m(rb' 


r  ) 

a 


=  A„  .  f.  .  (r,  , 

£s  j  £s  j  ,  m  b 


(4.3-15) 


where  the  spectroscopic  coefficient  A  .  includes  such  quantities 

36  S  ]] 

as  fractional  parentage  coefficients  and  the  interaction  strength. 

Summarizing  the  results  to  this  point,  eqn.  (4.3-1)  can  now 
be  written  as 
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where  the 

1  /  2 

(2 j+1) 


"reduced"  amplitude,  3  ,  is  given  by; 

i^  mamh  =  (-1) sb  mb <&sm,m  -  m  I j  M  -M  > 

sj  a  b  1  J  B  A 


and  where, 


x  < 


S  s  m  -m J  s  m  -  in  >(2j+l) 
a  b  a  b  a  b 


l  / 2 

Ps  j 


(4.3-17) 


(2  j+1)  1/2  iVm 

sj 


/<3?Jd?bxJ(_)  (£b4>f 


.  (r  ,r  )  x  ^  (k  ,r  ) 
£s] , m  baAa  a  b 


(4.3-18) 


Referring  to  eqn.  (4.3-18)  we  see  that  the  integration  over 
dr^  and  dr^  is  a  six  dimensional  integral,  which  in  general  is  a 
difficult  task.  To  overcome  this  difficulty  the  zero-range 
approximation  is  introduced.  In  this  approximation  particle,  b, 
is  assumed  to  be  produced  at  the  same  point  that,  a,  is  absorbed. 
Thus 


r 
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(4.3-19) 
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(4.3-20) 


The  resultant  integration  is  then  reduced  to  a  three  dimensional 
integral . 

It  has  been  shown  by  Austern  et.  al_.  ,  (Au  64)  that  the  zero- 


range  approximation  is  valid  if  only  small  momentum  components 
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of  the  distorted  wave  are  important,  characteristic  of  a 
surface  interaction. 

If  we  now  incorporate  the  zero-range  approximation  and 
the  partial  wave  expansion  of  eqn.  (4.3-8)  into  the  calculation 
of  eqn.  (4.3-18)  we  get,  upon  realizing  the  orthogonality  of 
the  spherical  harmonies: 


where 


8im  (6)  =  l  P1"  p“  (0)  fT”sj  ;  m>o 

3  L  L  L  Lk  Lb  LaLb 

a  b  a  b 
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(4.3-21) 


and 


<  L  £00  L  oxL,  £  m-m  L  o> 
b  1  a  b  1  a 


(4.3-22) 


ny-  M  M 

f IjS L  =  FT-r  /xL  (kb-rr)f£sj(r)  XL  <ka'r>dr  (4.3-23) 

ab  abA  B  B  J  a 

and  Pm(0)  is  the  associated  Legendre  polynomial  of  order  L,m. 

Two  principal  features  characteristic  of  direct  reactions 

are  realized  in  eqn.  (4.3-21),  viz: 

(i)  The  shape  of  the  angular  distribution  is 

given  by  an  infinite  sum  of  Legendre 

£  s  "i  . 

polynomials.  However  f  ,  quickly  tends 

L  L 
a  b 

to  zero  for  L  >  k  R  and  L,  >  kR  (where 

a  a  b  b 


R  is  the  nuclear  radius) .  The  resultant 


truncated  series  then  establishes  the 
characteristic  stripping  shape. 
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(ii)  The  parity  of  the  transition  is  given  by 
£ 

(-1)  as  dictated  by  the  requirement 

that  L  +L  +  £  be  even  in  the  Clebsch- 
a  b 

Gordon  coefficient,  <L  l  oo ! L  o>. 

b  8 

The  differential  cross-section  for  unpolarized  projectiles 
and  unpolarized  target  nuclei  is  given  by 


da  _  ya  pb  kb  Z It |  2 _ 

( 2 Trh 2 )  2  ka  (2J  +1)  (2S  +1) 

A  a 

where  and  y^  are  the  reduced  masses  of  the  respective  pairs, 
and  the  sum  is  over  M^,  m^,  Mg ,  and  m^ .  Upon  substituting  for 
T(eqn.  4.3-16);  summing  over  M  ,  m  ,  M  ,  and  rn  ;  and  realizing 

A  8  id  JD 

the  orthogonality  properties  of  the  Clebsch-Fordon  coefficients, 
eqn.  (4.3-24)  becomes: 
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(4.3-25) 


It  will  be  noted  that  in  the  resultant  expression,  eqn.  (4.3-25), 
the  contributions  from  £,  s,  and  j  are  all  incoherent.  This  is 
not  the  case  when  the  effects  of  spin-orbit  coupling  are  consi¬ 
dered.  Then  only  the  j  summation  is  incoherent  (Sa  64) . 

The  results  developed  to  this  point  are  applicable  to 
direct  reactions  in  general  (i.e.  stripping,  pickup,  inelastic 
scattering,  etc.).  Let  us  now  specialize  the  treatment  to  the 
single  nucleon  stripping  reactions  characterized  by,  A(d,n)B, 
i.e.: 
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a 

(4.3-26) 

'V 

-> 

(4.3-27) 

rb 

->  r 

n 

V 

+  V  . 

(4.3-28) 

np 

The  choice  of  interaction  potential,  V  ,  (i.e.  the 

np 

neutron-proton  interaction)  is  based  on  the  following  approxi¬ 
mation  . 

As  viewed  in  the  recoil  channel,  the  total  interaction 
potential  is  given  by: 


Vm  =  V  r  =  V  +  V  . 
T  nB  np  pA 


(4.3-29) 


The  residual  interaction,  V,  is  given  by  subtracting  off  the 
optical  potential,  U^B : 


V=V  +  V  -  U  ; 
np  pA  . nB 


(4.3-30) 


and  is  approximated  by  the  central-force  potential 


V  =  V 


np 


(4.3-31) 


The  nuclear  matrix  element  of  eqn.  (4.3-1)  can  now  be 


written  as 


J  M 
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np  A  A  d  np  p  n 

where  ^  represents  the  internal  coordinates  of  A;  a  the  spin 
A  P 

struecture  of  the  captured  proton;  and  a  the  spin  structure  of 
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the  emitted  neutron.  The  Jacobian, 
transformation  from  the  coordinates 
and  is  given  by 


J,  corresponds  to  the 

r  and  r  to  r  and 
pA  np  n 


J  = 


M,  M 
d  a 


M  (M.+M  ) 
p  d  A 


(4.3-33) 


To  evaluate  the  matrix  element  of  eqn .  (4.3-32)  the 

following  forms  for  the  wave  functions  are  usually  taken. 


(i)  The  nuclear  wave  function  $  is  expanded  into  a 
fractional  parentage  expansion  over  the  eigenstates  of  the 
target  nucleus  and  the  single  particule  state  of  the  captured 
proton ,  i .e  .  : 
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(4.3-34) 


The  function  may  be  regarded  as  proportional  to  the  single 

particle  wave  function  of  the  captured  proton.  It  can  be 
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separated  into  its  radial,  orbital,  and  spin  components  in  the 
following  way, 


^BA’(? 


PA' 
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(4.3-35) 


The  number  n  equals  the  number  of  equivalent  nucleons  in  the 
orbital  and  is  introduced  to  antisymmetrize  the  wave  function. 
The  coefficient,  a0  .  ,  is  the  fractional  parentage  coefficient 
and  the  product, 


S0  .  =  n(a0  .)  2  (4.3-36) 

Zj  Zj 

is  the  spectroscopic  factor.  It  is  a  measure  of  the  probability 
that  the  final  nucleus  is  made  up  of  the  configuration  of  the 
target  nucleus  and  the  captured  particle  in  the  single  particle 
state . 

(ii)  The  neutron  is  considered  as  an  elementary  particle 
with  no  internal  structure,  and  is  given  by  the  normalized  spin 
function : 

<t>%r%(a).  (4.3-37) 

(iii)  Only  the  S  state  structure  of  the  deuteron  is 
considered,  which  makes  it  factorable  into  the  product  of  a 
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Hulthen  wave  function  and  a  spin  function,  i.e.: 


1  mi  ( r  ,  o  ,  o  )  =  N 
np  p  n 


■  -ar  -fir 

e  np  .  e  p  np 
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(4.3-38) 


where  the  normalization  coefficient  is  given  by, 


N  = 


gg ( a  +  3) 

tt  (  3  -  g) 


(4.3-39) 


The  parameter,  g,  is  given  by: 

/MeT 

g  =  - 9_  (4.3-40) 

li2 

where  M  equals  the  average  of  the  neutron  and  proton  mass  and 
e  is  the  binding  energy  of  the  deuteron  (e  =  2.2245  Me  V.) . 

The  parameter  3  is  usually  taken  as 


3  =  7  a. 


(4.3-41) 


Using  the  above  forms  for  the  respective  wave  functions 
the  matrix  element,  eqn .  (4.3-32)  can  now  be  evaluated.  Since 

the  interaction  potential,  does  not  depend  on  the  co¬ 

ordinates,  r  A,  the  orthogonality  condition  demands  that  the 

only  contributing  state  of  the  target  nucleus  is,  J'  =  J  and 
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M'  =  M  .  Thus,  upon  expanding  the  spin  structure  of  the  deuteron 
into  the  spin  structure  of  the  neutron  and  proton  and  coupling 
the  composite  spin  to  the  value,  s,  eqn.  (4.3-32)  can  be  written 
as : 

£  S  -m 

J<B ,n I V  | A, d>  =  J  E  i'  (-1)  n<JAjMa,MB-MA|jBMB> 

£s  j 

X  <«.sm,md-mn  I  jMB-MA><sdsnmd-mn  I  smd-mn> 


*  2S+1  n%jP£j(rPA> 


(4.3-42) 


where  m  =  M„  -  MA  +  m .  -  in . 

B  A  d  n 


Upon  comparing  eqn.  (4.3-42)  with  eqn.  (4.3-14)  we  see 


that , 


/2  S  n  +  l 


m 


n  a0  .  p0.(r  A)Y*  ( r  )  V(  r  )  <J>  *  (r__)  (4.3-43) 


r£s  j  ,  m  2S  + 1  £j  £jv“pAy  £  np'  np,Tdv  np 


where  consistent  with  the  ideas  of  eqn.  (4.3-15)  we  make  the 
separation : 


A 


/2  Sd+1  ^ 

2S+1  n  V^rnp^d  ^rnp^ 


(4.3-44) 


( r  )  . 

pA 


(4.3-45) 


The  zero-range  approximation  is  introduced  by  writing. 
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(4.3-46) 


with 


D  =  |  V  (r)  c p  (r) dr. 
o  *  np  a 


(4.3-47) 


are  usually  taken  as  eigensolutions  of  the  Schrddinger  equation 

for  a  potential  well  of  the  Wood-Saxon  type  (eqn.  (4.3-4)). 

The  parameters  used  (V  ,  r  ,  and  a)  are  usually  those  taken  from 

o  o 

proton  elastic  scattering  work. 

In  the  event  that  the  stripped  proton  is  captured  into  a 
nuclear  orbital  of  defined  spin  and  parity  charge,  the  summa¬ 
tion  over  £  and  s,  in  eqn.  (4.3-39),  is  superfluous  (i.e., 
s  =  ^  and  £  =  j  -  h  as  dictated  by  the  parity  change)  .  Also 
from  the  structure  of  eqn.  (4.3-42)  we  see  that  the  following 
selection  rules  apply: 


£ 

(i)  the  parity  change  is  given  by  (-1)  as  dictated  by 
the  presence  of  the  Y ™ ,  and 


(ii)  the  final  spin  value,  J  ,  as  limited  by  the  Clebsch- 

B 

Gordon  coefficients,  is: 


J  -£-  ^ I  <  J  <  J  +£+^ . 
A  “  B  a 


(4.3-48) 


(ii)  The  D.W.B.A.  Angular  Distributions 


A  program,  written  by  B .  E.  F.  Macefield  (Ma  64)  and 
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embodying  the  theory  of  the  preceeding  section,  was  used  to 

calculate  the  angular  distributions  corresponding  to  the 

formation  of  the  ground,  0.335,  1.751,  and  2.433  levels  of 
2  1 

Na  via  the  Ne20(d,  n)Na21  reaction.  A  local,  Wood-Saxon 
optical  potential  (eqn.  4.3-3)  (with  zero  volume  absorption 
and  spin-orbit  terms)  was  used  to  evaluate  the  scattering  state 
wave  functions.  In  addition  the  zero-range  approximation  was 
used  throughout. 

The  average  optical  potential  parameters  given  by  Perey 
and  Buck  (Ro  66)  were  used  for  the  neutron  wave  functions. 
However,  the  corresponding  parameters  for  the  deuteron  were  not 
as  easily  obtained.  A  literature  search  revealed  that  no 
optical  model  analysis  had  ever  been  done  for  the  elastic 
scattering  of  deuterons  off  Ne2 0  and  that  the  closest  relevant 
work  was  some  elastic  scattering  data  at  10.95  MeV.  (Ta  60) 
and  11.6  MeV. (Ja  65).  Working  under  the  assumption  that  the 
energy  difference  between  the  elastic  scattering  data  and 
the  reaction  data  could  be  accounted  for  by  a  small  variation 
in  V,  the  above  data  was  subsequently  analysed  by  an  optical- 
model  search  code,  written  by  C.  M.  Perey  (0b66).  A  total  of 
five  optimum  fits,  ranging  in  V  value  from  58.0  MeV.  to  170.0 
MeV. ,  were  obtained  for  the  two  sets  of  data.  (The  correspond¬ 
ing  variation  in  the  other  parameters  was  not  significant.) 

When  the  aforementioned  sets  of  deuteron  parameters  were 
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substituted  into  the  D.W.B.A.  program  little  change  in  shape 
was  predicted,  however,  a  20%  change  in  magnitude  was  noted. 

Since  this  produced  a  corresponding  20%  variation  in  the 
spectroscopic  factors  it  became  necessary  to  single  out  a 
unique  set  of  deuteron  parameters.  This  was  achieved  by 
appealing  to  the  hypothesis  that  the  real  part  of  the  optical 
potential,  for  a  deuteron,  should  be  near  the  sum  of  the  neutron 
and  proton  potentials.  This  would  place  it  near  100  MeV. 

Deuteron  parameters  of  near  equal  value,  for  the  two 
different  sets  of  data  (10.95  and  11.6  MeV  data),  were  then 
averaged  and  the  results  closest  to  a  real  potential  value  of 
100  MeV  was  chosen.  The  results  so  selected  corresponded  to  a 
real  potential  value  of  121.995  MeV. 

In  an  effort  to  improve  the  D.W.B.A.  fits  a  search  through 
V  was  made.  This  was  consistent  with  the  idea  that  the  difference 
in  energy  between  the  scattered  and  reaction  data  could  be 
accounted  for  by  a  variation  in  V.  A  marginal  degree  of  im- 
provement  was  achieved  by  using  a  real  potential  value  of  95.0 
MeV. 

Figs.  4-4  and  4-5  illustrate  respectively  the  results  of 
the  elastic  scattering  analysis  and  the  D.W.B.A.  calculations, 
while  Tables  4-6  and  4-7  summerize  the  potentials  used. 
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Table  4-6 

Optical-Model  Parameters  used  by  Scattered  Waves 


Parameter 

Deuteron 

Neutron 

V 

121.995;  95.00 

MeV. 

48.0  -  0.29E 

MeV. 

W 

D 

25.721 

MeV. 

9.6 

MeV. 

R 

r 

1.228 

f . 

1.27 

f . 

R. 

1 

1.291 

f . 

1.27 

f . 

a 

r 

0.671 

f . 

0.66 

f . 

a . 

i 

0.625 

f . 

0.47 

f . 

R 

c 

1.300 

f . 

- 

Table  4-7 


Optical-Model  Parameters  for  Captured  Proton 


Parameter 

gnd . 

0.335 

1.751 

2.433 

V 

56.22 

55.59 

65.14 

50.28 

R 

1.25 

1.25 

1.70 

1.25 

r 

a 

0.65 

0.65 

0.65 

0.65 

r 

R 

1.30 

1.30 

1.30 

1.30 

c 
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Referring  to  Fig.  4-5  we  see  that  only  the  0.335;  2.433; 
and,  to  a  lesser  extent,  the  1.751  levels  are  characterized  by 
the  stripping  mechanism  (and  thus  validate  the  corresponding 
spin  and  parity  assignments)  whereas,  the  ground  state  is  not 
so  described.  A  valid  determination  of  the  shape  of  the  1.751 
level  is  hampered  by  a  very  low  cross  section  (0.563  mb/st.  at 
95.86°) . 

Spectroscopic  factors  for  the  0.335,  1.751,  and  2.433 
levels  are  given  in  Table  4-8  and  were  obtained  from  the  equation, 

H  (0)  =  1-48  sfTI  S  0(6,s.p.  (4.3-49) 

2J  +1 

where  a(0)s  p  is  the  output  value  given  by  Macefield's 

program.  The  factor  1.48  comes  from  the  Hulthen  wave  function. 


Table  4-8 

Spectroscopic  Factors 


Deuteron 

Potential 

Nuclear 

Level 

Max . 

da  ,  , .  x 

dfi  (exp  t') 

mb, 

vst' 

Max . 

m  (Ma-> 

S 

0.335 

6.03 

(43.88°) 

16.21 

(44.0°) 

0.37 

V,  =  121.995 
d 

1.751 

0.563 

(95.86°) 

4.90 

(85.0°) 

0.11 

2 . 433 

50.05 

(0°) 

198.32 

(0°) 

0.25 
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Table  4-8  (continued) 


Deuteron 

Potential 

Nuclear 

Level 

Max. 

da  ,  , ,  v 

dQ  (eXP  <=•> 

<— ) 

Max. 

do  (Ma-> 
<§> 

S 

0.335 

6.03 

(43.88°) 

15.82 

(42.0°) 

0.38 

V,  =  95.00 
d 

1.751 

0.563 

(95.86°) 

4.69 

(90.0°) 

0.12 

2.433 

50.05 

(0°) 

209.59 

(0°) 

0.24 

4 . 4  Discussion 

(i)  Approximations  made  in  the  D.W.B.A.  analysis 

Referring  to  Fig.  4-5,  we  see  that  even  for  those  levels 
exhibiting  stripping  structure  the  fits  are  far  from  ideal. 

There  is  justification  for  this  inconsistency^  however,  when 
one  realizes  the  many  approximations  made  in  the  D.W.B.A.  theory, 
viz  : 

(i)  no  contribution  due  to  compound 

reaction  formation  was  considered; 

(ii)  only  first  order  perturbation  theory 
was  used  (eqn.  4.3-1); 

(iii)  no  volume  absorption  or  spin-orbit 
effects  were  included  in  evaluating 
the  scattering  state  wave  functions; 

(iv)  the  optical  potentials  giving  rise  to 
the  scattering  state  wave  functions 
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were  assumed  to  be  local  potentials; 

(v)  the  zero-range  approximation  was  used 
throughout ; 

(vi)  the  interaction  potential,  V,  was 

reduced  to  a  central-force,  neutron- 
proton  interaction; 

(vii)  only  the  S  state  structure  of  the  deuteron 
was  considered;  and 

(viii)  the  captured  proton  radical  wave  functions, 
u£j (r  )  r  were  taken  as  solutions  of  a 
spherically  symmetric  potential  character¬ 
istic  of  shell  structure,  whereas,  in  the 
case  of  deformed  muclei,  such  as  Na2 1 ,  a 
non-spherical  potential  is  required. 

Of  the  approximations  listed,  number  (i)  can  contribute 
the  greatest  error  to  both  the  shape  and  magnitude  of  the 
angular  distribution.  For  in  reality,  the  true  transition 
amplitude  consists  of  two  terms,  viz: 

T  =  Td  +  Tc  (4.4-1) 

where  T  refers  to  direct  reaction  formation  and  T  to  compound 
reaction  formation.  Then  referring  to  eqn.  (4.3-24),  the 
differential  cross  section  becomes 
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da 


Z 

i 


cos 


d).  T1T1 
DC 


(4.4-2) 


(where  i  refers  to  the  spin  substates) .  Thus  if  compound 
reaction  formation  is  not  negligible,  the  effects  of  lTc|2 
and  the  interference  terms  can  alter  significantly  the  shape  and 
magnitude  of  the  angular  distribution. 

Point  number  (viii)  is  the  second  source  of  major  error, 
and  results  in  a  spectroscopic  factor  smaller  than  would  be 
expected  from  the  Nilsson  wave  function  formalism.  The  remaining 
approximations  tend  to  be  of  secondary  effect. 

( i i )  The  absence  of  stripping  structure  in  the 
formation  of  the  ground  state  of  Na2 1 

The  most  identifying  feature  of  the  results  of  Fig.  4-5 
is  the  absence  of  stripping  structure  in  the  formation  of  the 
ground  state  of  Na2 1 .  On  the  basis  of  the  D.W.B.A.  theory  of 
section  4.3  the  shape  of  an  angular  distribution,  for  single 
nucleon  capture  into  a  state  of  well  defined  spin  and  parity, 
is  determined  primarily  by  the  Z  value  of  the  captured  nucleon. 
Since  the  ground  state  of  Na2 1  has  a  well  defined  spin  and  parity 
value,  characteristic  stripping  structure  should  be  evident. 

One  reconciles  this  anomaly  by  appealing  to  the  following 
arguement . 

If  the  state  excited  is  poorly  defined  as  a  single  particle 
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level  a  small  spectroscopic  factor,  and  thus  direct  reaction 
component  of  the  transition  amplitude,  is  predicted.  The 
presence  of  any  compound  reaction  formation  then  becomes  sig¬ 
nificant  and  assumption  (i) ,  of  the  previous  section  breaks 
down.  The  generalized  expression  for  the  differential  cross 
section,  eqn.  (4.4-2),  must  then  be  employed,  with  a  resultant 
variation  in  the  shape  of  the  angular  distribution. 

One  can  say  qualitatively  from  the  poor  agreement  between 
the  theoretical  and  experimental  angular  distributions,  for 
the  ground  state,  that  the  compound  reaction  component  of  the 
transition  amplitude  appears  to  be  the  dominate  term.  This 
arguement  may  well  apply  to  the  1.751  level  also. 

( i i i )  Collective  model  interpretation  of  Na21 

The  collective  model  interpretation  has  been  applied  with 
a  high  degree  of  success  to  Na2 1  and  other  nuclei  in  the  ld-2s 
shells:  (Bh62) ,  (Ch63) ,  and  (Ho65) .  In  the  case  of  NaZi  the 

deformed  core  consists  of  10  protons  and  10  neutrons  coupled  to 
an  equilibrium  shape  of  a  prolate  spheroid.  The  remaining 
single  proton  is  thought  of  as  orbiting  this  core  under  the 
influence  of  a  non-spherical  potential. 

The  Hamiltonian  describing  such  a  system  is  given  by, 


=  H 


+  H 


H 


rot 
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(4.4-3) 
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where  H  refers  to  the  collective  rotational  motion  of  the 
rot 

deformed  core  and  refers  to  the  single-particle  motion  of 
the  extra-core  nucleons.  In  explicit  form  these  terms  are 
given  by: 


H 

rot 


J  ( J+l) 


21 3 


(K-ft) 


2 


and 


21 


+  J 


H 

o 


(Z 

P 


V 


(4.4-4) 


(4.4-5) 


where  J  is  the  total  spin  of  the  nucleus,  j  is  the  composite 
spin  of  the  extra-core  nucleons,  K  is  the  projection  of  J  onto 
the  body  fixed  3-axis,  ft  is  the  projection  of  j  onto  the  body 
fixed  3-axis,  and  I  is  the  moment  of  inertia  of  the  deformed 
core.  The  rotational-particle  coupling  (R.P.C.)  operator, 

J+_.  +  J__.  +  ,  serves  to  couple  together  states  of  different  K 

value,  namely  AK  =  ±1  (AK  =  0  if  K  =  %) .  The  sum  in  eqn.  (4.4-5) 
is  usually  reduced  to  just  one  particle  by  dropping  the  inter¬ 
particle  coupling  term  ?p*3q* 

If  one  ignores  the  R.P.C.  term  of  eqn.  (4.4-4)  and  makes 
the  assumption  of  a  rotationally  symmetric  deformation,  i.e., 

K  =  ft  (a  feature  characteristic  of  low  lying  levels)  one  gets 
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for  the  resultant  properly  normalized  and  symmeterized*  solution 
of  the  composite  system. 


where 


JM^7Ta> 


2  J+l 


16tt  2 


J 


^7ra> 


+  (-1)J  ^  „  |  -fiiT a> 

M-iZ  ' 


H  ^7ra>  =  fiiTa> . 
o'  il 1 


(4.4-6) 


(4.4-7) 


The  eigensolutions ,  [fiTra>,  form  what  are  known  as  the 
Nilsson  model  wave  functions,  where  tt  is  the  parity  of  the  single 
particle  level  and  a  is  the  Nilsson  orbital  number.  In  the  shell 
model  limit  the  potential  used  by  Hq  becomes  spherically  sym¬ 
metric  and  the  corresponding  eigensolutions  become  the  degenerate 
shell  model  wave  functions  parameterized  by  y  (the  strength  of 
the  l2  splitting  term)  and  k  (the  strength  of  the  spin-orbit 
coupling) .  For  the  Nilsson  wave  functions  the  additional 
parameter,  ri ,  is  introduced  as  a  measure  of  the  deformation;  and 
since  a  well  defined  axis  exists  for  deformed  nuclei,  namely 
the  symmetry  axis,  the  spin  degeneracy  is  removed  and  n(=  K) 
becomes  a  constant  of  motion. 

In  evaluating  the  wave  functions  |  ^7ra>  a  spherical-limit, 

*  One  symmeterizes  the  wave  function  with  respect  to 
reflection  in  the  equatorial  plane. 
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shell  model  expansion  is  usually  made,  i.e.: 


|ft7Ta>  =  E  <N£j  |  ftTra>  |N£j>  (4.4-8) 

n£  j 

where  |N£j>  form  a  basis  set  of  shell  model  wave  functions 
characterized  by  the  phonon  number  N,  orbital  angular  momentum 
£,  and  composite  spin  j;  and  where  <N£j|ftTra>  are  the  expansion 
coefficients.  The  summation  in  eqn.  (4.4-8)  is  usually  restricted 
to  a  single  oscillator  shell,  i.e.,  one  N  value. 

With  the  above  forms  taken  as  eigensolutions ,  the  eigen¬ 
values  of  eqn.  (4.4-3)  become, 


E  ( J  ,  ft )  —  + 


21 


J(J+1)  -  2ft2  + 


6^a(-1) 


J  +^> 


(  J+%  ) 


(4.4-9) 


where,  a,  known  as  the  decoupling  parameter,  is  given  by 


a  =  -  £  <N  £  j  |  ftTT  a  >  2  ( j  +h )  (-1)^+2.  (4.4-10) 

j 

The  values  have  been  calculated  numerically  and  are  sum¬ 
marized,  for  orbits  of  interest,  by  the  Nilsson  diagram  of  Fig. 
4-6  (Ho6  5 )  . 

From  eqn.  (4.4-9)  we  see  that,  for  a  given  ft(ft^s),  the 
rotational  bands  are  built  upon  the  respective  single  particles 
levels  with  spin  predictions  given  by  J  =  ft ,  ft  +  1 ,  ft  +  2 ,  ... 

With  this  knowledge  the  following  interpretation  of  the  level 
structure  of  Na2 1  can  be  made. 
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Fig.  4-6:  Nilsson  Diagram 
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A.  J .  Howard  et  al. ,  (Ho  65)  have  estimated  the  defor¬ 
mation  parameter  for  A=21  nuclei  as  being,  n  =  4.0.  With 
reference  to  Fig.  4-6  we  see  that,  with  four  nucleons  going 
into  each  level,  the  Na21  nucleus  can  then  be  described  by  a 
single  proton  in  the  |  ftTra>  =  1 3/2  +  7>  state.  This  gives 
the  required  ground  state  spin  value  of  3/2+.  The  5/2+  (0.335) 
and  7/2  (1.751)  levels  can  then  be  interpreted  as  the  ft+1 
and  ft+2  elements  of  the  rotational  band  built  upon  the  j  3 / 2 7  > 
state.  The  l/2+  (2.433)  level  is  subsequently  interpreted  as 
a  particle  excitation  to  the  | 1/2  +  9>  state. 

To  this  point  the  analysis  has  neglected  the  R.P.C.  term 

•  tt  -f-  + 

which  would  then  mix  together  components  of  the  ft  =3/2  ,  1/2  , 

+ 

and  5/2  states.  However,  it  has  been  shown  by  Kelson  and 

TT  + 

Levinson  (Ke  64)  that  the  ft  =3/2  state  predominates  in  the 
3/2+,  5/2+,  and  7/2+  levels  of  Na2 1  to  the  amount  of  0.98, 

0.85,  and  0.80.  The  neglect  of  this  term  is  then  justified. 

In  the  context  of  the  D.  W.  B.  A.  theory  of  Section  4.3, 
a  spectroscopic  factor  represents  the  probability  to  which  the 
final  state  wave  function  can  be  interpreted  as  consisting  of 
an  extreme  single  particle  level  coupled  to  an  unexcited, 
spherically-symmetric ,  target-nucleus  core.  Such  an  inter¬ 
pretation  is  only  valid  for  nuclei  near  closed  shells.  For 
deformed  nuclei,  such  as  Na2 1 ,  one  must  account  for  possible 
rotational  and  vibrational  motion  of  the  core  in  addition  to 
introducing  non-spherical  potentials  for  the  extra-core  nucleons. 
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Following  the  outline  of  Satchler  (Sa  58)  the  theory  is 
modified  in  the  following  way. 

Referring  to  eqn.  (4.3-32)  we  see  that  under  the  assump¬ 
tion  that  v=v(rnp)  only,  the  quantity  of  interest  becomes  the 
overlap  integral  between  the  target  and  residual  nuclear  wave 
functions ,  i . e . : 

^«<rpA'?)  M  (5)  d?  =  .1  Vm(rpA)<JAjMAralJBV 

B  B  ^  A  A  jZ  JJ  ^ 

(4.4-11) 


where  once  again  we  have  made  a  fractional  parentage  expansion 
and  representing  the  wave  function  of  the  captured  proton. 
And  as  before  the  spectroscopic  factor  is  given  by, 


S0  .  =  n3?  .  .  (4.4-12) 

Zj  Zj 

In  evaluating  the  overlap  integral  of  eqn.  (4.4-11)  we  make 
the  following  substitutions  for  the  nuclear  wave  functions: 


M  (?)  = 

jama 


~2J  +1 

h\ 

A 

_16tt2 

J  -j  J 

I  .  .  /  -i  \  A  A  A 

dm  k  |Q7KX>a  +  (-1}  dm  -K 

A  A  A  A 


-^7Ta> 


1 


^Vib 


( 4 . 4-13a) 


y  (£  r'  )  = 
jbMb  pA 


2JB+l" 

h 

^ISit2 

_ i 

D  B  Iftirar'  > 
mbkb  pA  B 


+  (-1)  B  3b  DmB  |-nirar;A>B 

B  B  r 


Cb- 


( 4 . 4-1 3b) 


(For  the  sake  of  generality  we  have  assumed  that  K  ^  S2,  and 
have  also  accommodated  any  vibration  structure  via  4>vib-) 
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Before  preforming  the  integration  of  eqn.  (4.3-11)  one  must 


realize  that  states  j Prar  >  are  related  to  the  body  fixed 

1  pA  B 

coordinate  system,  whereas,  the  function  b . _  (r  )  is  related 

j  Jim  pA 

to  the  space  fixed  system.  This  inconsistency  is  brought  into 


line  by  realizing  that  * 


£m  pA 


I 

m 


D J  . 
mm 


6 


Jim 


(r  ) 
pA 


Thus  the  integral  of  interest  is, 


(4.4-14) 


I  J  MJ  *  [-£ 

m'  x  J  M  pA 
B  B  ^ 


H  <5) 
A  A 


DJ  M 

mm  ' 


(4.4-15) 


Satchler  has  shown  (Sa  58)  that  the  results  of  such  an 


integration  yield. 


/S.  (K  fi  KB)  =  g 
J  £  A  A  B  B 


2  J  +1 
A 

2  J  +1 
B 


1J5 


<J  J  +  K  ,  ft  ~  ft  ,j  J  K  > 
A  A  B  A  1  B  B 


x  <  N  £  j  j  Q  -  ft^iraxc^l  $A>  (4.4-],r 

where  g  =  /2  if  either  K  =  ft  =  0  or  K  =  ^  =  0,  and  unity 

otherwise.  The  -  is  taken  according  to  whether  the  captured 
proton  enters  an  orbit  with  ft  =  |  ftg  -  |  .  The  selection  rule 

orovided  by  the  Clebsch— Gordon  coefficient  is; 

K  -  K  =  Q  -  !2d.  (4.4-17) 

A  B  A  B 

Since  most  ground  states  are  characterized  by  =  ftA 


r- 


)  (  t  3  ■  i  /}.  •  .  • 


' 


’  'r-  :c 


-72- 


this  requires  that  Thus  stripping  to  highly  excited 

states,  characterized  by,  K  /  ,  would  consequently  be 

retarded . 

In  the  situation  corresponding  to  the  Ne 2 c (d , n) Na2 1 
reaction  a  proton  of  defined  j , i  is  captured  by  an  even-even 
nucleus  for  which  =  I<^  =  £2  =0.  Eqn .  (4.4-16)  then  reduces 

to , 

SjJt  =  I  <  N£j|nTra>2  (4.4-18) 

where  the  overlap  integral  between  the  vibration  states, 

<(J)  |({)  :>,  has  been  set  equal  to  unity  consistent  with  the  inter- 

pretation  of  the  low-lying  structure  of  Na2 1 

Thus  by  exciting  successive  levels  of  a  given  rotational 
band  one  is  able  to  explore  the  corresponding  expansion  coeffi¬ 
cients  of  the  Nilsson  model  wave  function.  And  since 
<N£j|Q7Ta>2  represents  the  probability,  that  in  a  measurement, 
the  eigenstate  |J2ira>  will  be  found  in  the  state  |N£j>,  the 
classical  interpretation  of  a  spectroscopic  factor  still 
applies. 

Upon  applying  eqn.  (4.4-18)  to  the  spectroscopic  data 
of  Table  4-8,  one  gets  (realizing  that  the  ld-2s  shells  are 
characterized  by  N=2) : 


<2 

2 

5/2  | 

3/2  +  7 > 2  = 

1.12 

(4.4-19) 

<2 

4 

7/2  j 

j  3/2  +  7>  2  = 

0.46 

(4 . 4-20) 

<2 

0 

1/2  | 

1/2+9  > 2  = 

0.25 

(4.4-21) 
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An  important  conclusion,  relating  to  the  Ne 2 0  (d ,n) Na 2 1 
experiment,  can  now  be  made. 

Since  the  basis  vectors  |N£j>  form  a  complete  set  of 
functions  for  a  given  N  value, 

l  <N£j  |ftTra>2  =  1.  (4.4-22) 

The  results  of  eqns .  (4.4-19)  and  (4.4-20)  then  indicate  that, 

within  the  framework  of  a  purely  direct-reaction  formulation, 
the  measured  cross  sections  for  the  0.335  and  1.751  levels  are 
too  high.  This  anomaly  can  be  attributed  to  the  lack  of 
attention  paid  to  the  compound  reaction  component  of  the  transi¬ 
tion  amplitude. 

If  one  makes  a  number  of  simplifying  assumptions,  the 
results  of  eqns.  (4.4-19),  (4.4-20),  and  (4.4-22)  can  be  used 

to  estimate  a  minimum  value  for  the  average  compound-reaction 
contribution  to  the  angular  distributions  of  the  0.335  and  1.751 
levels . 

First,  let  us  assume  that  the  direct  reaction  and  compound 
reaction  components  of  the  transition  amplitude  add  incoherently. 
Second,  let  us  assume  that  the  |ftTra>  =  1 3/2+7>  state  is  made  up 
of  only  the  |N£j>  =  \2  2  5/2>  and  \2  4  7/2>  basis  vectors. 

(This  corresponds  to  assuming  that  the  Na21  ground  state  is 
formed  wholly  via  the  compound  reaction  mechanism  and  that 
higher  £  value  components  of  the  basis  set  are  negligible.) 

Our  expansion  of  interest  is  then: 
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3/2+7 >  =  <2  2  5/2  |  3/2  +  7>  |  2  2  5/2> 

+  <2  4  7/2  |  3/2+7 >  |  2  4  7/2>,  (4.4-23) 


for  which 


<22  5/2  |  3/2  +  7 >2  +  <24  7/2  |  3/2  +  7>2  =  1.  (4.4-24) 

If  we  now  assume  that  the  average  compound-reaction  component 
is  a  constant  value,  for  both  the  0.335  and  1.751  levels,  we 
get  (using  the  results  of  Table  4-8,  eqn.  (4.4-18)  and  eqn. 
(4.4-24) )  : 


2 ( 5/2 ) +1 

6 . 03-c” 

,  2  ( 7/2 ) +1 

0 . 563-c 

2 

16.01 

+  2 

i 

• 

for  which  the  average  background,  compound-reaction  component 
is : 

c  -  0.6  mb/st  (4.4-26) 

This  result  then  forces  one  to  interpret  the  1.751  level 
as  being  formed  entirely  by  the  compound  reaction  mechanism, 
and  also  attributes  a  10%  compound  reaction  contribution  to 
the  0.335  level. 

Of  the  assumptions  made,  in  this  rather  over-simplified 
analysis,  that  of  an  equivalent  compound  reaction  contribution 
to  both  the  0.335  and  1.751  levels  is  perhaps,  in  first  thought, 
the  least  intuitively  justified.  One  might  argue  that,  since 
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the  £  value  for  the  captured  proton  is  4  for  the  1.751  level 
P 

and  only  2  for  the  0.335  state,  momentum  transfer  probabili¬ 
ties  would  be  greater  for  the  0.335  level  than  for  the  1.751. 
Such  an  argument  is  consistent  with  the  idea  that  the  proton- 
target  interaction  distance,  r  ,  must  be  greater  for  £^=4 

capture  than  for  £  =2  (as  dictated  by  the  classical  condition 

P 

that  £  =k  r  ) . 


P  P  'pA 


Such  an  argument  is  indeed  valid,  within  the  framework 
of  the  direct  reaction  formalism,  and  serves  to  account  for 
the  reduced  cross-section  of  the  1.751  level.  However,  within 
the  compound  reaction  formalism,  we  are  involved  with  the 
complete  absorption  of  the  deuteron  followed  by  the  subsequent 
emission  of  a  neutron.  That  is,  a  new  point  of  view  is  taken 
and  the  angular  momentum  arguments  must  be  altered  accordingly. 
The  conditions  of  interest  here  are  the  conservation  of 
channel  spin  and  parity. 

If  we  assume  that  the  compound  nucleus,  Na2 2 ,  is  formed 
primarily  by  S-wave  aeuterons  (which  is  consistent  with  the 
idea  of  the  complete  absorption  of  a  rather  large  nucleus  at 
a  bombarding  energy  of  3.924  MeV)  then  the  entrance  channel 
spin  and  parity  value  are: 


(4.4-27) 
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The  conservation  of  channel  spin  allows,  respectively  for  the 
ground,  0.335,  1.751,  and  2.433  levels,  angular  momentum  values 
of : 


L 

n 

(gnd)  = 

i — 1 

o 

,  2,  3 

(4.4-29) 

L 

n 

(0.335) 

=  1, 

2,  3,  4 

(4.4-30) 

L 

n 

(1.751) 

=  2, 

3,  4,  5 

(4.4-31) 

L 

(2.433) 

=  1, 

2 

(4.4-32) 

for  the  neutron  exit  channel,  as  dictated  by: 

|J  21  -  S  -  1  I  <  L  <  J  2i+Sn+l  (4.4-33) 
1  Na  n  1  -  n  -  Na^ 

However,  since  all  four  levels  are  positive  parity  states  the 
conservation  of  channel  parity, 

nf  =  (-l)Ln  nn  1^2  i.  (4.4-34) 

restricts  the  L  values  to: 


L^(gnd)  =  0,2 

(4.4-35) 

L  (0.335)  =  2,4 
n 

(4.4-36) 

L  (1.751)  =  2,4 
n 

(4.4-37) 

L  (2.433)  =  2. 

(4.4-38) 

With  these  results  we  see  that  the  presence  of  any  centri- 
f ugual-barrier  effects  must  influence  the  0.335  and  1.751  levels 


equally.  Thus  the  assumption  of  an  equivalent  compound  reaction 
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contribution,  to  the  0.335  and  1.751  levels,  is  given  foun¬ 
dation.  In  addition  to  this  conclusion,  if  one  couples  the 
results  of  eqn.  (4.4-35)  and  the  realization  that  the  stripping 
cross  section  is  proportional  to,  2J  +1,  (eqn.  4.3-25)  the 
following  argument  can  be  fabricated  in  support  of  the  assumption 
that  the  ground  state  is  formed  wholly  via  the  compound  reaction 
mechanism. 

With  a  ground  state  spin  value  of,  Jf  =  3/2,  we  see  that 
the,  2Jf+l,  term  would  tend  to  retard  the  direct  reaction  for¬ 
mation  of  this  state  with  respect  to  the  0.335  level.  (Both 

levels  have  £  =2  values,  therefore,  angular-momentum  transfer 

P 

arguments  apply  equally.)  However,  the  allowability  of  a 

L  (gnd)  value  of  zero  would  tend  to  enhance  the  compound  reaction 
n 

formation  of  this  state  over  that  attributed  to  the  0.335  level. 
The  complementary  effect  of  these  two  restriction  tends  to 
validate  the  assumption. 

Let  us  now  turn  our  attention  to  the  2.433  level.  Having 
excited  only  the  J=Q  level  of  the  \h+9>  band  it  is  difficult 
to  form  any  conclusions  about  the  appl.icibility  of  the, 

l  <n£ j  |  *2+9 > 2  =  1,  (4.4-39) 

£  j 

rule.  However,  realizing  that  highly  excited  states  are 

generally  characterized  by  K  /ft ,  the  selection  rule  of 

eqn. (4.4-17)  would  demand  small  spectroscopic  factors  for  the 


' 

. 

(,  '•  .o  r  ‘  1 


....  '  '  '  ' 


78- 


higher  energy  members  of  the  \h+9>  band.  The  value  of, 

<20  \  |  h+9>2  =  0.25,  for  the  2.433  level  is  then  probably 

quite  realistic.  This  result  then  implies  a  strong  direct 

reaction  contribution  to  the  formation  of  this  level,  which  is 

further  substantiated  on  the  basis  of  angular-momentum  transfer 

arguments.  With  a  L  (2.433)  value  of  2,  the  compound  reaction 

n  1 

contribution  is  probably  no  greater  than  that  attributed  to  the 
0.335  and  1.751  levels.  This  would  represent  a  1.1%  contribution. 

One  surprising  feature  about  the  2.433  level  is  that,  in 
spite  of  the  fact  that  it  corresponds  to  a  single  particle  exci¬ 
tation,  it  is  poorly  represented  by  a  spherical ly-svmmetric 
shell  model  wave  function  (as  indicated  by  a  spectroscopic 
factor  of  only  0.25).  It  may  be  that  the  K=ft  selection  rule  is 
coming  into  play  at  this  energy. 
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APPENDIX  A 


ASSOCIATED  PARTICLE  ANALYSIS:  ASPANL 
A-l  Introduction 

ASPANL  is  a  combined  FORTRAN  IV  and  MAP*  language  program  written  for 
an  IBM  7040  computer.  Its  primary  function  is  to  compute  the  ratio,  R, 
between  the  number  of  neutrons  incident  on  the  neutron  detector  to  the 
number  of  neutrons  produced  in  coincidence  with  a  collimated  flux  of  the 
kinematically  associated  particles.  Fig.  A-l  is  a  schematic  outline  of 
the  experimental  configuration. 

The  following  quantities  are  required  in  the  determination  of  R: 

(i)  the  size,  shape,  and  spacing  of  the  charged  particle  collimators 
C2  and  C3, 

(ii)  the  size  and  nature  of  the  target:  T, 

(iii)  the  size  and  shape  of  the  incident  beam  collimator:  Cl, 

(iv)  the  distance  of  the  neutron  detector,  D2 ,  from  the  target,  and 

(v)  the  kinematic  inter-relationship  of  the  reaction  particles. 

These  variables  interact  in  such  a  way  as  to  produce  a  2-dimensional 
distribution  for  the  associated  neutron  flux.  The  quantity,  R,  is  then 
determined  by  computing  the  ratio  of  the  sum  of  this  distribution  over  the 

Macro  Assembly  Program 
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FIGURE  A-l  Associated  Particle  Geometry 
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neutron  detector  to  the  total  sum  of  the  distribution.  The  basic  task  of 


ASPANL  is  then  to  determine  the  distribution  of  the  associated  neutron  flux. 

This  is  achieved  by  the  following  synthesis  of  the  experimental 
situation.  First  one  defines  a  matrix  to  represent  the  distribution  and, 
in  terms  of  the  mathematics  of  the  program,  positions  it  at  the  same 
distance  and  angle  as  the  neutron  detector.  The  second  initialization 
feature  is  to  divide  the  "active  volume"  of  the  target  into  a  specified 
number  of  differential  volume  elements.  The  program  then  computes  the 
coordinates  for  each  volume  element  and  asks  the  question:  for  a  given 
element  is  it  possible  to  pass  a  straight  line  through  collimators  C2  and 
C3?  If  so,  the  range  of  angular  variations,  in  both  the  vertical  and 
horizontal  planes,  are  noted  and  the  corresponding  kinematically  associated 
variations  in  the  neutron  angle  are  computed.  A  count  is  then  added  to 
those  elements  of  the  matrix  subtended  by  the  associated  neutron  angular 
variations.  Each  element  is  interrogated  in  this  way  and  the  resultant 
contents  of  the  matrix  represent  the  2-dimensional  flux  distribution. 

To  minimize  computational  time  the  more  redundant  features  of  this  inter¬ 
rogation  cycle  are  written  in  MAP  machine  language. 

The  approximations  made  in  the  above  synthesis  of  the  problem  are: 

(i)  A  constant  reaction  cross  section  is  assumed  across  the 
angular  variation  allowed  by  the  collimators. 

A  quadratic  approximation  is  made  to  the  kinematic  inter¬ 
relationship  of  the  two  particles  over  the  angular  variation 
allowed  by  the  collimators . 


(ii) 
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(iii)  The  number  of  matrix  elements  subtended  by  two  different  volume 
elements,  which  in  turn  subtend  equal  solid  angles  by  the 
collimators,  may  differ  due  to  the  different  angular  orienta¬ 
tions  of  the  volume  elements  and  the  non-linear  kinematic 
relationship  between  the  two  particles.  Thus,  as  a  consequence 
of  adding  a  unit  count  to  each  matrix  element  subtended,  the 
contribution  to  the  neutron  distribution  will  be  greater  for 
that  volume  element  which  subtends  the  greater  number  of 
matrix  elements.  Whereas,  in  reality,  each  volume  element 
should  contribute  in  a  manner  proportional  to  the  solid  angle 
subtended  by  the  collimators. 

These  approximations  are  employed  in  an  effort  to  minimize  computational 
time  and  thus  any  small  error  that  might  be  introduced  is  vindicated.  In 
the  case  of  point  (iii)  the  point  of  view  adopted  by  the  program  results 
in  a  significant  saving  in  machine  time,  for  one  is  then  able  to  program 
the  more  redundant  features  of  the  interrogation  cycle  in  fixed  point. 

The  program  is  written  with  sufficient  versatility  to  accommodate  a 
number  of  different  situations,  viz,  one  can  either  specify  the  distance 
of  the  neutron  detector  or  have  it  computed  by  an  algorithm  which  positions 

it  at  such  a  distance  that  it  just  begins  to  subtend  the  constant  flux 

% 

portion  of  the  distribution.  In  addition  either  a  circular  or  rectangular 
aperture  in  collimator  Cl  can  be  specified.  However,  one  is  restricted  to 
circular  apertures  only  in  collimators  C2  and  C3.  Additional  restrictions 
are  that  the  target  be  of  finite  depth  and  that  the  face  of  the  neutron 
detector  be  of  circular  geometry. 
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The  primary  application  of  the  program  has  been  to  a  circular  gas 
target  cell  and  the  D(d,n)He3  reaction.  However,  its  application  is 
easily  extended  to  other  reactions  provided  a  similar  target  geometry  is 
employed . 

The  above  serves  as  a  general  introduction  to  the  operation  of  ASPANL. 
For  a  more  detailed  description  refer  to  (Bu  67). 

A- 2  Data  Input 

All  input  information  to  ASPANL  is  entered  via  punched  cards.  Each 
card  contains  a  two  digit  control  code  and  up  to  eight  pieces  of  data. 

Each  control  code  directs  the  program  in  compliance  with  the  respective 
input  information. 

The  following  table  describes  the  contents  of  each  card.  The  format 
used  throughout  is: 

FORMAT (IX, 12 , IX, 8E9 . 0) 


TABLE  A-l  INPUT  DATA 


Code 

DATA 

01 

SI,  S2 ,  Rl,  R2 ,  XSIZE,  YSIZE,  ZSIZE,  DIAM 

51  =  distance  from  the  center  of  the  target  cell 

to  the  first  collimator,  C2  (inches) 

52  =  distance  from  the  center  of  the  target  cell 

to  the  second  collimator,  C3  (inches) 
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Table  A-l  continued 


CODE 


DATA 


R1 

R2 

XSIZE 

YSIZE 

ZSIZE 

DIAM 


aperture  radius  of  the  first  collimator,  C2 
(inches) 

aperture  radius  of  the  second  collimator,  C3 
(inches) 

length  of  target  cell  (inches) 

width  of  collimator  Cl  (inches) 

height  of  collimator  Cl  (inches) 

diameter  of  the  neutron  detector  (inches) 

— variable  DIAM  is  optional  and  need  not  be 
included. 


Control  card  01  is  to  be  entered  if  a  rectangular  aperture  in  colli¬ 
mator  Cl  is  used. 


14 


SI, 


S2 ,  Rl,  R2 ,  XSIZE,  YSIZE,  DIAM 


SI 

““ 

distance  from  the  center  of  the  target  cell 
to  the  first  collimator  C2  (inches) 

S2 

= 

distance  from  the  center  of  the  target  cell 
to  the  second  collimator  C3  (inches) 

Rl 

= 

aperture  radius  of  the  first  collimator  C2 
(inches) 

R2 

= 

aperture  radius  of  the  second  collimator  C3 
(inches) 

XSIZE 

= 

length  of  the  target  cell  (inches) 

YSIZE 

= 

aperture  diameter  of  collimator  Cl  (inches) 

DIAM 

= 

diameter  of  the  neutron  detector  (inches) 

— variable  DIAM  is  optional  and  need  not  be 
included . 

Control  card  14  is  to  be  entered  in  exchange  for  control  card  01  if 
a  circular  aperture  in  collimator  Cl  is  used. 
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Table  A-l  cont inued 


CODE 

DATA 

02 

Ml,  M2,  M4 , 

Q,  Tl,  LARGE 

Ml 

mass  of  the  incident  particle  (A.M.U.) 

M2 

mass  of  the  target  particle  (A.M.U.) 

M4 

mass  of  the  primary  particle  (neutron)  (A.M.U.) 

Q 

Q  value  of  the  reaction  (MeV) 

T1 

energy  of  the  incident  particle  (MeV) 

LARGE  = 

a  variable  specifying  that  in  the  case  of  a 
double  valued  reaction  the  higher  (or  lower) 
energy  of  particle  3  (He')  is  used  in  evaluat¬ 
ing  the  neutron  flux  distribution,  i.e. 

LARGE  >_  -0.001  :  higher  energy  is  used. 

LARGE  <  -0.001  :  lower  energy  is  used. 

LARGE  unspecified  :  higher  energy  is  used. 

Control 

card  02  need  not  be  given  if  control  card  10  is  given. 

03 

PRES,  TEMP, 

CURRNT 

PRES 

pressure  of  gas  within  the  target  cell  (cm  Hg) 

TEMP 

temperature  of  gas  within  the  target  cell  (°C) 

CURRNT  = 

the  beam  current  incident  on  the  target  cell 
(yamps) 

If  this  card  is  entered,  the  target  thickness  is  computed  assuming  a 
deuterium  target. 

See 

also 

control  cards 

05  and  08. 

04 

NUMPTS 

NUMPTS  = 

the  number  of  differential  volume  elements 
into  which  the  "active  volume"  is  to  be 
divided.  From  this  number  the  incremental 
distances  DELX,  DELY,  DEYZ  are  computed. 

NUMPTS  is  set  to  500.0 

if  neither  control  card  04  or  07  is  specified. 
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.Table  A-l  continued 


CODE 

DATA 

05 

THETAO ,  SIGMA 

THETAO  = 

the  central  angular  position  of  the  collima¬ 
tor  system  C2  and  C3  with  respect  to  the  beam 
direction  (degrees) 

5°  £  THETAO  175° 

SIGMA 

the  differential  cross  section  of  the  reac¬ 
tion  corresponding  to  particle  3  (He3) 
being  located  at  THETAO  (mbarns/sterd. ) 

If  control  card  03  is  entered,  a  reaction  count  rate  is  computed 
assuming  the  D(d,n)He3  reaction. 

If  control  card  03  is  not  entered,  SIGMA  is  ignored  and  the  parameter 
UNTCNT  is  set  equal  to  unity. 


(BLANK) 

Cont rol  card  06  causes  the  first  78  characters  of  the  next 
card  to  be  accepted  as  title  information. 


07 


DELX,  DELY 

DELX 


DELZ 


=  incremental  distance  along  length  of  active 
volume  (inches) 


DELY  =  incremental  distance  along  width  of  active 
volume  (inches) 


DELZ  =  incremental  distance  along  height  of  active 
volume  (inches) 


This  card  overrides 


08 


UNTCNT 


UNTCNT  = 


control  card  04. 


the  count  rate  per  cubic 
within  the  target 


inch  per  steradian 


This  card  overrides  the  value  of  UNTCNT  calculated  from  the  data 
entered  via  control  cards  03  and  05. 


*  0AT3HT 

/.MO  1 2 

b>>;  j\.  U03  Q Jt  X'  fon  roi  lOBf  T  6  t  '■  3  -  fa!5i  ?>  iO'j 

.no  7  B’31  9H(n«b)Q  srf3  snimuaes 


.ylj  u  o.i  Inupa  3aa  ai  PfcOTHU 

— 


(3MAJSI 


5o  nal  '  sDatielb  ipnx  - 

(barfonP)  STiuXov 


(feaxionl)  antflov 

.AO  biBD  loi lnoo  zabiusvo  birBD  airfT 


Table  A-l  continued 


CODE 


DATA 


09 


DN 


DN  =  the  distance  from  the  center  of  the  target 

cell  to  the  neutron  detector  (meters) 

If  DN  is  not  entered,  it  is  computed  subject  to  the  condition  that 
the  neutron  detector  just  begin  to  subtend  a  constant  flux. 

1.016  <_  ON  £  6.0  (m) 

If  DIAM  is  not  given  or  if  all  six  angles  are  given,  via  control  card 
10,  then  DN  is  set  to  the  longest  value  (6.0  m)  if  it  is  not  speci¬ 
fied. 


10 


THETAP,  THETAO ,  THETAM,  ZETAP,  ZETAO ,  ZETAM 


THETAP  = 


THETAO  = 


THETAM  = 


ZETAP 


ZETAO 


ZETAM 


minimum  angle  allowed  by  collimator  C2  and 
C3  (degrees) 

central  angle  of  collimators  C2  and  C3 
(degrees) 

maximum  angle  allowed  by  collimators  C2 
and  C3  (degrees) 

the  neutron  angle  kinematically  associated 
with  THETAP  (degrees) 

the  neutron  angle  kinematically  associated 
with  THETAO  (degrees) 

the  neutron  angle  kinematically  associated 
with  THETAM  (degrees) 


All  these  angles,  except  THETAO,  are  normally  calculated  from  the 
data  entered  via  control  card  01  or  14. 

The  data  entered  via  control  card  10  overrides  all  previously  calcu¬ 
lated  values. 


11 


NORM 


NORM 


a  variable  whose  magnitude  determines  whether 
the  matrix  representing  the  neutron  flux  dis¬ 
tribution  is  normalized  with  respect  to  count 
rate  on  the  basis  of  the  computed  target 
thickness,  cross  section,  and  beam  current. 
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Table  A-l  continued 


CODE 

DATA 

NORM  >  -0.001  :  the  scatter  diagram  is  normalized 

NORM  unspecified  :  the  scatter  diagram  is  normalized 

NORM  <_  -0.001  :  the  scatter  diagram  is  not  normalized 

12 

(BLANK) 

This  control  code  terminates  data  input  and  initiates  execution  of 

the 

run. 

13 

(BLANK) 

This  control  code  terminates  the  program  by  calling 


SUBROUTINE  EXIT 


The  data  cards  are  to  be  entered  in  order  of  increasing  code  number 
(except  for  code  14  which  replaces  01)  and  all  parameters s  except  those 
on  control  cards  07,  09,  and  11  remain  unchanged  from  run  to  run.  Table 
A-2  is  a  sample  listing  of  input  data  illustrating  the  need  for  only  a 
portion  of  the  available  input  cards. 

The  three  cards  following  control  card  13  form  the  calling  sequence 
for  AUTOPLOT,  a  library  plotting  program  written  by  J.  F.  Easton  (Ea  64). 
Its  relationship  to  ASPANL  is  discussed  under  heading  A-3  Data  Output. 
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TABLE  A- 2  SAMPLE  INPUT  DATA 

$ ENTRY 

14  0.60 

4.60  0.038  0.100  1.00  0.182  3.50 

02  2.014740 

2.014740  1.008986  3.266  2.96 

03  3.50 

20.0  1.5 

04  8000.0 

05  31.2 

3.50 

06 

D(D,N) HE- 3 

E  2.96  THETA  31.2  MARCH  16  1967  RUN  NO.  8 

09  2.012 

12 

13 

$ EXE CUTE 

AUTOPT 

PLOT  ALL 

PLOT  DONE 

^  >v  [ 


T40TUA 


A- 11 


A- 3  Data  Output 

The  output  of  each  run  consists  of  a  data  sheet,  a  scatter  diagram, 
and  a  graphical  plot,  represented  respectively  by  Figs.  A-2,  A-3,  and  A-4. 

The  first  portion  of  the  data  sheet  is  a  kinematic  analysis  of  the 
reaction.  The  quantities  Ml,  M2,  M3,  and  M4  are  respectively  the  masses 
of:  the  incident  particle,  the  target  particle,  the  associated  particle 

(He3),  and  the  primary  particle  (n)  with  Tl,  T3,  and  T4  their  respective 
energies.  The  three  angles  listed  are  respectively:  the  minimum,  central, 
and  maximum  angles  allowed  by  collimators  C2  and  C3.  In  the  event  of  a 
double  valued  reaction  the  corresponding  angles  and  energies  of  the  lower 
energy  neutron  group  are  printed  alongside  the  angles  and  energies  of  the 
higher  energy  group. 

The  next  portion  under  the  heading  "Given  Data"  is  just  a  restatement 
of  some  of  the  input  data. 

The  quantities  under  the  heading  "Computed  Data"  includes:  a  state¬ 
ment  of  the  number  of  volume  elements  employed,  based  on  the  incremental 
values  of  DELX,  DELY,  and  DELZ;  a  statement  of  the  time  taken  to  compute 
the  flux  distribution;  and  an  evaluation  of  the  quantity  UNTCNT,  i.e.  the 
count  rate  per  cubic  inch  per  steradian  within  the  target. 

The  final  portion  under  the  heading  "Scatter  Diagram  Data"  contains 
the  important  results  of  the  program.  The  quantities  XNDEL  and  YNDEL  are 
respectively  the  width  and  height  of  the  matrix  elements  making  up  the 
flux  distribution.  These  quantities,  when  multiplied  by  the  dimensions 
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of  the  matrix  (101x101) ,  give  the  total  width  and  height  of  the  scatter 
diagram.  The  next  sequence  of  numbers  relates  to  the  flux  distribution 
and  outlines,  among  other  things,  the  total  count  rate  (counts/sec) 
represented  by  the  distribution.  The  quantity  entitled  "Integrated 
Count  in  the  Primary  Detector"  is  the  count  rate  incident  on  the  detector, 
and  the  quantity  "Ratio  of  Total  Counts  to  Those  in  the  Detector"  is  the 
all  important  parameter,  R,  discussed  earlier.  The  remaining  data  is 
self  explanatory  and  serves  to  compare  the  solid  angle  of  the  neutron 
detector  to  the  solid  angle  of  the  flux  distribution. 

The  next  output  sheet  (Fig.  A-3)  is  a  scatter  diagram  of  one-half  of 
the  2-dimensional  flux  distribution.  Only  one-half  of  the  distribution  is 
shown  since  it  is  symmetrical  about  the  bottom  line.  The  statement  "Matrix 
Was  Normalized"  refers  to  the  fact  that  the  distribution  was  normalized  to 
an  upper  limit  of  10,  represented  by  the  stars.  Lower  valued  integers 
then  correspond  to  a  decrease  in  intensity.  With  a  little  imagination  it 
is  possible  to  visualize  a  3-dimensional  image  from  this  figure. 

The  remaining  output  sheet  (Fig.  A-4)  is  a  graphical  plot  of  the  flux 
distribution  obtained  by  a  call  to  AUTOPLOT.  The  x-axis  measures  the  radius 
of  the  distribution  and  the  y-axis  the  count  rate  per  square  inch.  The 
dotted  curve  is  a  plot  through  the  center  of  the  distribution  in  the  hori¬ 
zontal  plane.  (It  corresponds  to  the  bottom  row  of  Fig.  A-3.)  The  curve 
made  up  of  crosses  is  a  plot  through  the  center  of  the  distribution  in  the 
vertical  plane.  The  two  bars,  located  symmetrically  about  the  zero  radius 
position,  correspond  to  the  placement  of  the  neutron  detector  and  are 
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determined  internally  in  accordance  with  the  value  of  DN,  the  distance  of 
the  neutron  detector  from  the  target.  It  will  be  noted  that  the  title 
entered  via  control  card  06  appears  on  each  data  sheet. 
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